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NORMAS PARA ADMISSÃO 


DE TRABALHOS É 


INSTRUÇÕES PARA OS AUTORES 


1. A Revista TÉCNICA, em princípio, publicará anualmente 
quatro números de natureza técnico-cientifica. 


2. Atendendo ao espaço disponivel e carácter científico que se 
deseja imprimir à Revista, as três formas de publicação previstas, 
revestir-se-ão das seguintes características: 


ARTIGOS POR CONVITE — Trabalho de síntese, no dominio 
da investigação científica ou no domínio do desenvolvimento, 
contemplando actividade desenvolvida entre nós, ou pelo menos, 
de interesse para o desenvolvimento científico e tecnológico do 
Pais. 


ARTIGOS TÉCNICO-CIENTÍFICOS — Trabalho de investiga- 
ção científica ou de aplicação tecnológica a um domínio bem 
definido da ciência ou da técnica e de que o autor seja um especia- 
lista. 

NOTAS CIENTÍFICAS — Destinam-se essencialmente à apre- 
sentação de novos resultados no dominio da investigação cienti- 
fica ou do desenvolvimento tecnológico. Pretende-se uma nota 
curta, realçando os resultados alcançados e referenciando conve- 
nientemente trabalhos afins, evitando a exposição longa de ante- 
cedentes, que deverá ser remetida para referências. 

3. Os artigos por convite e os artigos técnico-científicos, deverão 
ser escritos em língua portuguesa. As notas científicas poderão ser 
escritas em língua estrangeira, de preferência em inglês. 


4. Os autores deverão atender na elaboração dos artigos às se- 
guintes condições: 

— os artigos convidados não poderão exceder 30 páginas dacti- 
lografadas, incluindo o espaço reservado a ilustrações; 

— os artigos técnico-científicos não poderão exceder 15 páginas 
dactilografadas, incluindo o espaço reservado a ilustrações; 

— as notas científicas não poderão exceder 8 páginas dactilogra- 
fadas. 

5. Os trabalhos devem ser dactilografados a dois espaços, em 
folhas modelo Ad, devidamente numeradas. 


6. Os trabalhos devem ser enviados em duplicado, recomendan- 
do-se aos autores que conservem uma cópia, 
7. Conjuntamente com o trabalho, devem ser enviados; 


— título do trabalho em português e inglês; 

— resumo € abstract; 

— nome, direcção e telefone do autor e demais dados que facili- 
tem o processo de revisão de provas. 


8. Os autores deverão assegurar-se, antes de enviarem os origi- 
nais, que estes não contenham erros, falta de dados, etc., pois que, 
na fase de revisão de provas, apenas se corrigirão erros de compo- 
sição tipográfica. 


9. As referências bibliográficas serão indicadas no texto por or- 
dem numérica, usando algarismos árabes entre colchetes, e a lista 
das mesmas, no final do artigo, obedecerá às seguintes normas: 


Artigos de revistas: nome(es) do autor(es) e iniciais nome da re- 
vista, número do volume, página e ano, entre parêntesis. O título 
do artigo só deverá ser mencionado se isso constituir matéria 
fundamental para esclarecimento. 

Livros: nome do autor e iniciais, título, editora, ano de publica- 
ção e páginas, 


10. Os Quadros ou Tabelas, numerados em romano, devem figu- 
rar no texto datilografado, podendo ser encimados por um título 
conciso. 


11. Os gráficos devem ser desenhados a tinta da china, em folhas 
separadas e prontos para serem reproduzidos. Quaisquer letras 
ou algarismos que neles figurem, deverão ser suficientemente 
grandes para permitir a sua redução que, em princípio, fica ao 
critério do editor. 

As folhas com os gráficos, figuras ou fotografias, serão numera- 
das, indicando-se o trabalho a que dizem respeito. Além disso, no 
texto, indicar-se-à o local para a respectiva inclusão. 


NOTA: 4 TÉCNICA salienta a necessidade de os autores, 
quando for o caso, obterem em tempo as devidas autorizações de 
reprodução de excertos de outras obras, sendo esta matéria da 
inteira responsabilidade do autor. 


A TÉCNICA procurará, sempre que possível e por princípio, o 
parecer de individualidades de reconhecida competência técnico- 
científica, no sentido de apurar da originalidade, qualidade e inte- 
resse científico dos trabalhos que lhe são enviados. Esta atitude 
não se aplica aos artigos convidados, a menos que os autores o 
desejem. À Direcção da TÉCNICA, caberá sempre a ceitação ou 
recusa dos originais que lhe sejam enviados. A TÉCNICA notifi- 
cará os autores ou devolverá os trabalhos submetidos a aprecia- 
ção num prazo máximo de 90 dias após a recepção do original. Os 
trabalhos publicados, serão registados com a data de recepção do 
original na TÉCNICA e disso se dará pública notícia em respec- 
tiva nota de fim de página. 

A TÉCNICA oferecerá cinco exemplares da Revista, por cada 
trabalho publicado, independentemente do número de autores. 
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EDITORIAL 


As preocupações expressas nos nossos editoriais ainda se mantêm. 

A TÉCNICA continua a enfrentar limitações de ordem material que têm impe- 
dido a regularidade de publicação que se pretende e que, estamos, certos, se atin- 
girá. 

Este facto levou-nos a condensar numa única tiragem os 3.º e 4.º números de 
1984 e a meditar sobre métodos e formas a adoptar em futuras edições. 


No que respeita a original cumpre-nos registar a resposta positiva às nossas 
solicitações. Contudo, o que nos têm enviado é extenso e vai além das normas que a 
revista fixou. Sobre esta matéria não podemos deixar de salientar o preço por 
página exigido por algumas revistas estrangeiras. Por outro lado, e no que respeita 
a notas científicas nada se tem recebido. 


Tudo leva a crer que as dificuldades materiais sejam superadas mas o que se nos 
escapa é a existência de um número elevado e rico de notas científicas que é um dos 
aspectos mais importantes, não só no que respeita à projecção da revista mas 
também da imagem que dá da comunidade que serve. 


M. de ABREU FARO 


Homenagem ao Professor 
Ferreira de Macedo 


Estão em curso diligências no sentido de ser prestada homenagem ao Professor António Augusto 
FERREIRA DE MACEDO, que foi notável professor do Instituto Superior Técnico, de 1927 a 
1947, onde ensinou, magistralmente, Matemáticas Gerais e Geometria Descritiva, cadeiras básicas 
comuns a todos os cursos ministrados no IST. 

A Câmara Municipal de Lisboa dará todo o apoio às Homenagens, e está disposta a implantar um 
monumento em sua memória numa Praça de Lisboa, local que, quando o monumento lá estiver e 
for inaugurado, passará a chamar-se PRAÇA PROFESSOR FERREIRA DE MACEDO, 

Os fundos para a concretização de tal monumento serão obtidos por subscrição pública, para a 
qual contribuirão, espera-se, todos quantos, como alunos ou amigos, tiveram o privilégio de contac- 
tar ou conhecer tão excelsa figura de MESTRE e CIDADÃO. 

As eventuais contribuições monetárias, que antecipadamente se agradecem, devem ser concreli- 
zadas por transferência bancária para a “conta” HOMENAGEM AO PROF. FERREIRA E 
MACEDO, com o n.º 1685872/000/001/009, expressamente aberta no BANCO FONSECAS 
BURNAY, Dependência Urbana, Av. Fontes Pereira de Melo, 4, em Lisboa. 

Eventuais pedidos de esclarecimento ou sugestões, podem ser dirigidos para: 


“Comissão de Homenagem ao Prof. Ferreira de Macedo” 
Av. Estados Unidos da América, 82-3.º Esq. — 1700 LISBOA 


Professor Pedro Luis Borges Teixeira 
Por ocasião do primeiro aniversário 
do seu falecimento 


O Prof. Borges Teixeira nasceu no dia 5 de Janeiro de 1941 e faleceu em 
16 de Fevereiro de 1984. | 
Pelo seu valor intrínseco, pelo modo como generosamente pôs as suas altas 
qualidades ao serviço do Instituto Superior Técnico, com 
independência da tristeza e palavras de convenção que sempre 
se motivam numa vida curta, subitamente interrompida em pleno brilho, 
o Departamento de Engenharia Electrotécnica e de Computadores deliberou desde 
logo assinalar a obra do Prof. Borges Teixeira. 

Aliou-se a Técnica a esse propósito. 

Acontece que neste momento a Direcção da Revista recai sobre mim. 
Assim, em relação à justa homenagem que esta Revista presta 
ao Prof. Borges Teixeira, não posso deixar de acrescentar um pouco mais, 
consubstanciado no facto de ter sido seu professor, de termos colaborado, 
de nos ligar uma amizade muito grande e antiga. 
Sentimentos desta natureza não se constroem no abstracto ou em inexplicável 
emotividade mas radicam-se, isso sim, em afinidades de ideal universitário. 
Não me compete julgar no pormenor, disso se ocupou em excelente artigo 
a Secção de Electrônica do nosso departamento. 
Uma última nota científica do Prof. Borges Teixeira, também publicada 
neste número, dá exemplo e notícia do seu modo próprio de analisar uma ciência 
que especialmente cultivou: a Electrónica. 
A par de uma sensibilidade e dotes de artista, o Prof. Borges Teixeira 
possuía uma inteligência lúcida e profunda, uma capacidade rarissima de aliar 
a prática eficiente do Engenheiro com uma teoria que dominava 
bem e pairava alto, 
Passou de modo efémero pela Indústria e nesse trânsito aproveitou o melhor. 
Estudioso incansável, possuidor de uma cultura universal, um sorriso 
permanente de disponibilidade acrescido, por vezes, com uma ironia sem mal que 
ainda era inteligência e aviso consagrou-se totalmente ao Técnico, 
deu-se integralmente à Escola. 
Desligado de interesses materiais, no espírito mais puro de um universitário, 
viveu nesse ideal e receava por esse ideal. 


M. ABREU FARO 
Fev. 85 
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RESUMO 


Neste trabalho abordam-se aspectos teóricos da 
mieracção electrão-alvo relevantes para a compreensão 
das caracteristicas da espectrocospra de perda de ener- 
gia de electrões, nomeadamente no que se relaciona 
com a observação de transições electrónicas e vibracio- 
nais opticamente interditas. São referidos diferentes 
tipos de dispositivos experimentais descritos na litera- 
tura e apresentam-se alguns espectros servindo de 
exemplo das potencialidales da técnica. 


ABSTRACT 


A brief review is given of the status of energy 
electron loss spectroscopy. The article discusses the 
physics of electron impact concerning the vibrational 
and electronic molecular excitations. Special attention 
to optical forbiden transitins is grven. Different types 
of spectrometers are described. Some spectre are pre- 
sented which points out the potentialities of the tech- 
nique. 


| — INTRODUÇÃO 


Designa-se por espectrometria de electrões todo um 
conjunto de técnicas baseadas na análise da energia 
cinética de electrões emitidos por uma amostra. Con- 
soante o tipo de radiação primária, ela pode tomar 
varias designações. Delas destacamos a espectroscopia 
do fotoelectrão e a de perda de energia de electrões. 

Na espectroscopia do fotoelectrão, o alvo é subme- 
tido a um feixe de radiação electromagnética, em geral 


Original recebido para publicação em 29/9/84. 


raiosX ou radiação ultravioleta. Não nos referiremos a 
ela neste trabalho, pois, por si só, constituiria tema 
para um longo artigo de revisão. Encaminhamos o lei- 
tor para as referências [1-5] como exemplo de traba- 
lhos publicados susceptíveis de lhe dar uma ideia do 
seu desenvolvimento. 

A espectroscopia por impacte de electrões, cuja 
abordagem constitui o objectivo deste trabalho, baseia-se 
na medida da energia perdida por electrões monociné- 
ticos nos diferentes processos de colisão em que inter- 
vêm. 

Os primeiros trabalhos no domínio da difusão de 
electrões reportanise ao início do século [6]. No 
entanto, foi só nos últimos vinte anos que esta técnica 
sofreu uma evolução apreciável. Isto foi devido ao 
desenvolvimento de novos métodos de selecção em 
energia na produção de electrões monocinéticos, ao 
aumento da precisão dos detectores e analisadores de 
partículas com carga eléctrica e ao incremento de novos 
métodos computacionais de aquisição e tratamento de 
dados. O interesse crescente por esta espectroscopia 
está igualmente associado às perspectivas que oferece 
nc estudo de numerosos processos atómicos ou mole- 
culares, pouco acessíveis pela espectroscopia óptica. 
É o caso da observação de transições opticamente inter- 
ditas, estados acima do nível de vácuo, nomeadamente 
os plasmões e ainda estados de iões negativos. Para 
além do apoio que fornece à caracterização estrutural 
dos átomos e/ou moléculas, é possível assinalar vários 
exemplos de aplicação em problemas concretos, como 
no domínio da identificação de substâncias. 

Na figura 1 está representado esquematicamente um 
tipo de espectro de electrões rectrodifundidos por uma 
superfície sólida [7]. Nele podemos distinguir três 
regiões: A região 4 contém um pico de perda nula, 
que corresponde à detecção de electrões difundidos 
elasticamente (é o pico elástico); a zona B recolhe as 
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— ENERGIA — 


Fig. 1: Espectro típico de electrões rectrodifundidos por uma 
uma superfície sólida: A — electrões difundidos elasti- 
camente (pico elástico), B — Zona de perda de energia, 
C — electrões secundários lentos ou secundários verda- 
deiros. 


contribuições dos electrões que estiveram na origem de 
excitações electrónicas, vibracionais e rotacionais das 
moléculas do meio irradiado; finalmente a parte € 
refere-se aos electrões secundários lentos ou secundá- 
rios verdadeiros. Em certos casos, é possível evidenciar 
nesta última região uma estrutura atribuível à emissão 
de electrões, por autoionização ou por efeito Auger. 
A análise da energia de electrões com esta origem cons- 
titui, aliás, a base de outras duas formas de espectros- 
copia de electrões: Auger e de Autoionização [3,8]. 

Neste artigo vamos incidir particularmente a nossa 
atenção sobre a zona B do espectro. A revisão que 
faremos começará por abordar aspectos teóricos das 
colisões inelásticas electrão-molécula e que permitem 
evidenciar as possibilidades de estudo que a técnica 
oferece. Depois de sintetizarmos algumas formas de 
concretizar experimentalmente a espectroscopia de perda 
de energia, apresentaremos alguns espectros obtidos por 
vários autores. Finalizaremos o artigo com a indicação 
de aplicações possíveis da técnica no domínio da Quí- 
mica e da Física. 


2 — ASPECTOS TEÓRICOS DA ESPECTROCOPIA 
DE PERDA DE ENERGIA DE ELECTRÕES 


Um estudo teórico da espectroscopia por impacte 
de electrões, terá que abordar necessariamente, o pro- 
blema da difusão inelástica de electrões por moléculas 
do meio irradiado. Tal abordagem é usualmente reali- 
zada separando o caso em que a energia dos electrões 
incidentes é suficientemente elevada de forma a que 
a sua velocidade seja muito superior à dos electrões 
K do meio (v > Ze'/h), daquele em que esta desi- 
gualdade não se verifica. Na primeira situação é possí- 
vel considerar a colisão como rápida e não ressonante: 
o electrão primário actua como uma perturbação impul- 
siva do sistema alvo que responde como se fosse um 
sistema isolado. No segundo, o electrão incidente pode 
ser fortemente acoplado ao alvo e torna-se indiscernível 
dos electrões orbitais. 
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Neste parágrafo vamos portanto separar estes dois 
tipos de interacções (colisões rápidas e colisões lentas) 
com o objectivo de prever a natureza das transições 
electrónicas produzidas por cada uma delas. As exci- 
tações de estados vibracionais serão igualmente tra- 
tadas. 


2.1 — Excitações electrónicas por difusão inelástica 
2.1.1 — Colisões rápidas 


A teoria da difusão inelástica de electrões por áto- 
mos foi elaborada por Bethe [9] no domínio da vali- 
dade da aproximação de Born (v > Ze*/h). Ela foi 
retomada por Lassetre [10], Brongersma [11] e Ino- 
kuti [12]. Citamos igualmente os livros de Massey [13] 
e de Bonham [14] bem como o artigo de revisão de 
Celotta e Huebner [15]. 

De acordo com a teoria de Bethe, a secção eficaz de 
uma colisão levar um átomo do seu estado fundamental 
(lunção de onda Y,) a um nível discreto caracterizado 
por uma energia de excitação E, (função de onda W,) é 


dada por: 
do 2 Ki -4 2 
q = tda da) (1) 


Nesta expressão a, = (W/m e)”, em que m e e são res- 
pectivamente a massa e a carga do electrão e 2 é o 
ângulo sólido. O vector de transferência de momento 


é definido por 


! 
y 


— K, sp Kr (2) 


e K; = |K; e Kr = |Kf são respectivamente os vecto- 
tores de onda de momento inicial e final. A quanti- 
dade |£,| representa, por sua vez, o valor absoluto do 
elemento de matriz de posição entre as funções de onda 
do estado inicial e final: 


Z 
n=<n|SeplianD|W> 6) 


| 


sendo Z o número total de electrões da malécula e r; 
o vector de posição do electrão j do alvo. Tendo pre- 


Fig. 2: Representação Requemátnica definindo a orientação dos 


vectores de onda K, e K, e transferência de momento a. 


= 0 = — 


sente a orientação relativa dos vectores Ki, K e q 
esquematizada na figura 2, é possível chegar a 


q =Ki+Ki—-2kK Kcos é (4) 


Por analogia com a teoria da absorção da radiação 
electromagnética, Bethe introduziu o conceito de força 
de oscilador generalizada (f.o.g.) o qual obedece à 
regra da soma > f = Ze se define pela expressão 


/ 2 
fn = E (q ao)? (e) (5) 


Sendo R = 13,605826 eV e a = 5,2917715. 10“ m. 
Tendo em conta a equação 5, a expressão da secção 
eficaz diferencial toma a forma 


do. E -. Ke l ho 
do + K; | (q do) | En/R (6) 


A secção eficaz diferencial depende portanto, simulta- 
neamente, das características do alvo, através de f,/En, 
e da partícula incidente pelo factor K$/K; q”. 

Para o objectivo do presente trabalho, interessa-nos 
considerar aqui o caso particular das «colisões rápi- 
das» acompanhadas de uma transferência de momento 
pequena. Neste caso o ângulo de difusão é reduzido, 
o que permite substituir 4 por 


q = (K; — Ky* + (K,09 (7) 
ou, fazendo intervir em vez dos momentos, a energia 


cinética E;, do electrão incidente e a energia de exci- 
teção da molécula alvo E,, 


"2 En E 
[e"+(2%) o 


Por sua vez, sendo K; — K; muito pequeno (K;= Kr), 
a equação 6 passa e ter a forma 


de RR L 

— e À —— RB 

do i BR dá 
Substituindo 8 em 9 obtém-se 

d RN E Í 

E4 ã& — = (10) 

dO E E e” + (E/2 E; 


No limite quando O” tende para zero, a secção eficaz 
diferencial (do/dQ), é dada por 


de asa E 
(ao )=16 Rn E; (11) 


mostrando que a probabilidade de excitação de um 
nível de energia E, depende da razão f /E| em condi- 
ções de pequena transferência de momento e de ângulo 
O próximo da difusão rasante. 

Se a transição Y, “WY, for opticamente permitida, 
pode demonstrar-se que a f.o.g. se reduz à f.o. óptica. 
Assim, o desenvolvimento em série do operador que 
figura na expressão 3 conduz a 


q q a” Sm ash 
meio mo SC ia à) dia 
o (iq a)” 
= A DES ari (12) 
ic m! 


um BAR a ma = a 
sendo q um vector unitário com a direcção q. Substi- 
tuindo 12 em 5 verifica-se que 


Lim ta = 00" — E ( ] 3) 
R 


q=—+ 0 


Nas condições limite definidas, a expressão 11 passa a 
escrever-se em função da força de oscilador óptica fj: 


do No ds 
( dO ). maça 16 do R E: E, (14) 


Podemos pois concluir que as transições electrónicas, 
produzidas por electrões rápidos em colisões inelásticas 
de pequena transferência de momento, devem obedecer 
às regras de selecção ópticas. Diz-se então que nos 
encontramos nas condições de validade da «aproxima- 
ção óptica», 

Se a transição Y, “” Yo for opticamente proibida é 
necessário tomar em consideração, no desenvolvi- 
mento 12, termos de ordem superior e explicar, deste 
modo, valores não nulos da força de oscilador. 

Em face do exposto podemos portanto concluir 
que um espectro de perda de energia de electrões rápi- 


SECÇÃO EFICAY 


ENERGIA (e) 


Fig. 3: (a) Espectro de perda de energia de electrões de 300 eV 
em benzeno (b) Espectro de absorção no U. V. do ben- 
zeno (Reprodução de [27]). 


dos, obtido segundo ângulos de difusão pequenos, deve 
coincidir, em termos de probabilidades relativas das 
diferentes excitações moleculares, com o espectro de 
absorção óptico. É o que se pode constatar na fig. 3, 
onde se representa um espectro do benzeno obtido com 
um feixe de electrões de 300 eV e segundo um ângulo 
de difusão razante [16] e um espectro de absorção no 
U. Y. [47). 

As expressões aplicáveis a transições entre estados 
discretos são facilmente modificáveis para transições 
para um estado final contínuo, substituindo em 9 o por 
do/dE e f, por df/dE. 

A probabilidade de excitação produzida por um 
electrão rápido em meios condensados é usualmente 
expressa em termos da constante dieléctrica, £, do meio. 
Tal procedimento fundamenta-se no facto da variação 
da parte imaginária de e, em função da frequência de 
excitação, apresentar uma estrutura que corresponde a 
pontos críticos nas transições de banda. Assim, na 
expressão da secção eficaz de excitação referente a uma 
transferência de energia E, = ho, o produto (1/,) 
df/dE escreve-se na forma [18]: 


1º df m i 
ho d(no) = dm* e? hºw im (- 4 (15) 
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sendo Im (— 1/£) a função de perda de energia defi- 
nida por 


tm ( E) si a (16) 
) = q+8 


em que N é o número total de moléculas do meio irra- 
diado, por unidade de volume e £, e £; representam, 
respectivamente, a parte real e imaginária da constante 
dieléctrica do meio. A abundância relativa de um estado 
dotado de uma energia de excitação ho é portanto defi- 
nida pela «função de perda», pelo que uma vez conhe- 
cidos os valores de £; e e podemos recolher informa- 
ções sobre a natureza das excitações primárias. 

A determinação de £, e £> pode ser feita a partir de 
medidas de reflectância de fotões [19, 20]. É o caso 
dos valores calculados para o poliestireno e represen- 
tados, na parte a, da figura 4, em função da energia 
de excitação [19]. A função de perda obtida a partir 
dos dados anteriores e indicada na figura 4b pode ser 


E lev) 


Elev) 


O 5 1 15 20 25 30 as 
A E (av) 


Fig. 4: (a) Valores experimentais de e, e e, do poliestireno 
obtidos por reflectividade em função da energia dos 
fotões [19], (b) Função de perda do poliestireno [20], 
(c) Espectro de perda de energia de electrões de 20 KeV 
difundidos em poliestireno [21]. 


comparada com um espectro de perda de energia de 


electrões rápidos em poliestireno [21], indicado na 


parte c, da figura 4. Verifica-se deste modo que a fun- 


ção de perda indica bem a probabilidade de excitação. 


(10783 em?) 


SECÇÃO EFICAZ TOTAL 


ENERGIA (eV) 


Fig. 5: Secção eficaz total de excitação de níveis moleculares 
do benzeno calculada segundo o modelo de Ochkur- 
“Rudge (Reprodução de [27]). 


Os métodos ópticos de obter a constante dieléctrica 
encontram, no entanto, algumas ificuldades experimen- 
tais, relacionadas, por exemplo, com o facto da reflec- 
tividade ser muito sensível ao estado da superfície, 
especialmente para radiações de baixo comprimento 
de onda (u.v. longíquo). Vários autores, baseados na 
possibilidade do espectro de perda de energia pode ser 
descrito pela teoria dieléctrica, têm procedido de forma 
inversa, isto é, determinando constantes ópticas a par- 
tir da análise dos espectros. Exemplos de tratamentos 
deste tipo, utilizando a expressão da dispersão de Kra- 
mers-Kroning são considerados nas referências [22-26] 
para substâncias orgânicas e inorgânicas. 


2.1.2 — Difusão inelástica de electrões lentos 


Se a velocidade dos electrões incidentes for vizinha 
ou inferior à dos electrões orbitrais, a aproximação de 
Born deixa de aplicar-se no estudo da difusão inelástica. 
A solução do problema teórico torna-se então mais 
complicada, pois é necessário considerar fenómenos 
ignoráveis nas «colisões rápidas»: processos de troca 
de electrões, polarização dos átomos em moléculas do 
alvo, distorção por estes últimos da onda plana asso- 
ciada ao electrão... 

Uma primeira extensão da aproximação de Born 
é a de Born-Oppenheimer na qual a função de onda, 
Y (r;,r») que caracteriza o sistema de dois electrões 
(incidente e difusor) é a função produto antissimé- 
trico [27]: 


Ear 
+* (nm)=N [e “w (n)+ 


+ et y, (19)] n (17) 


Nesta expressão n é a função de spin simétrica ou 
antissimétrica apropriada, N uma constante de norma- 
lização, r; e r; definem as posições dos electrões do 
meio e incidente, Y; (r;) é a função de onda do electrão 
do alvo, enquanto que e “2 é a do electrão incidente. 
A amplitude de difusão inelástica é dada, neste caso, 
por 


g> (0) = — 


2m pa ds iqr ” 
4 nº !) id fa 
E R | 1 tj diidi Í fa (Kim — Key) 4º (ti) 


2 2 
|-É is | Y, (ti) dridra (18) 


rj T; 


em que o sinal + se refere à excitação de um estado 
singuleto e o sinal — à de um nível tripleto. O pri- 
meiro termo da expressão corresponde à amplitude de 
difusão de Born e o segundo à difusão de troca. Este 
último, que explica o facto das amplitudes das transi- 
ções opticamente proibidas não serem nulas, foi poste- 
riormente reformulado por Ochkur [28, 29] e por 
Rudge [30, 31]. 

O tratamento de Ochkur-Rudge serviu de base a 
Matsuzawa [52] e Christophorou [27] para determi- 
narem secções eficazes de excitação do benzeno consi- 
derando-o como um sistema de 6 electrões 7 e um cor 
de electrões c. Não apresentaremos neste artigo as 
expressões obtidas, dada a sua complexidade e as 
mesmas não serem necessárias para o nosso objectivo. 
Assim limitamo-nos a reproduzir na fig. 5 os valores 
valculados para os estados singuleto e tripleto do ben- 
zeno e apresentados na referência [27]. Const'- -» 
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Fig. 6: 


Fig. 7: 
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o. 120º 150º 
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/ 


os | 30º 60º 


Distribuição angular normalizada correspondente à 
excitação dos níveis *E,, (T;), *B,, (T,) e *B,, (Ty). 
Energia incidente (a) 30 (eV) (b) 300 eV (Reprodução 
de [32]). 


E a, 


e e e o O E q O e 


(b) 


(a) Geometria de difusão de um electrão (b) Orienta- 
ção do momento dipolar associado com os osciladores 
adsorvidos. (Reprodução de [36]). 


portanto, que a probabilidade das transições optica- 
mente proibidas pode ter valores apreciáveis no caso 
de excitação por electrões lentos (E, < 50 eV). Assi- 
nale-se que esta probabilidade, para além de depender 
da energia incidente é igualmente função do ângulo de 
difusão O, como se pode observar na fig. 6. As transi- 
ções para o nível T:> (mais prováveis entre os estados 
tripletos) são favorecidos para valores pequenos de O. 
Na tabela seguinte são indicadas probabilidades máxi- 
mas das transições “Bj “ 'Ag, “Ent 'A, e “Bu T'Ar 
para vários valores da energia do electrão primário, 
com a indicação do ângulo de difusão a que se referem. 

Estes resultados mostram a importância da espec- 
troscopia de perda de energia de electrões lentos na 
observação de processos moleculares não revelados pela 
espectroscopia óptica. 


TABELA 


Valores máximos de secções eficazes de excitação dos esta- 


dos ºB,,, *E,, e ºB,, do benzeno /32| 
Bo dt IA, E, - A, 'B,. + IA, 
E, (eV) o Pre Tmnax Onix EU max PR Umax 
10 1351 0,223 49 0,861 
20 74 782.107 
50 58  3,93.107 27 0184 59 935.105 
100 50  4,13.10 
300 17º 478107 8 218103 17 LO7AO- 


2.2 — Excitações vibracionais por difusão inelástica 


Na sua fase actual, a espectroscopia vibracional por 
difusão inelástica de electrões está sobretudo orientada 
para o estudo de moléculas ou átomos adsorvidos em 
superfícies metálicas. Os trabalhos teóricos prevêem 
para este caso dois tipos de mecanismos de interacção 
electrão-vibração: as interacções dipolares e as interec- 
ções ressonantes. No que se segue tomaremos como 
ponto de partida dois artigos de revisão de Bertolini 
[33] e Lopez [34]. 


2.2.1] — Interacções dipolares 


Para distâncias entre o electrão incidente e uma 
molécula adsorvida elevadas, o electrão fica essencial- 
mente submetido ao campo eléctrico criado pela distri- 
buição de cargas formada pela molécula e pela carga 
de blindagem no substracto. A interacção sofrida pelo 
electrão é então de natureza dipolar [35, 36]. 


A geometria da difusão inelástica de um electrão 
por acção de um centro difusor está esquematizada na 
= = 


fig. 7a. Dafinem-se assim os vectores K; e Kr, bem 
como os ângulos O, 0º e 9” que fixam as orientações 


dos vectores em relação à superfície e ao plano de 
incidência. Na parte b da mesma figura representa-se 


o vector momento dipolar U do dipolo associado às 
oscilações do meio adsorvido, situado a uma distân- 
cia Zo da superfície e com uma orientação definida 
pelos ângulos O, e Ps. 

O sistema electrão-molécula-substracto pode ser 
definido por um hamiltoneano que para além da ener- 
gia cinética do electrão difundido, do potencial devido 
a superfície metálica limpa e de um termo correspon- 
dente aos modos vibracionais moleculares, contém 
ainda o potencial, V, devido à interacção do elec- 
trão com os campos dipolares das moléculas adsorvi- 
das [36]: 


du. r 
V=-elg+st— ju) — 19 
O primeiro termo nesta expressão provém do dipolo 
permanente e contribui só para a difusão elástica, 
sendo por isso desprezado. O segundo é devido ao 


dipolo oscilante na direcção de to e nele intervêm, a 
distância r do electrão ao centro de massa da molé- 
cula, s, uma coordenada normal de vibração e (du /ds) 
o momento dipolar dinâmico associado a esta coorde- 
nada. 

Partindo da equação (19) é possível deduzir a 
seguinte expressão para a secção eficaz total de difusão 
inelástica [36]: 


Em é [de 
E M pela 


2 
) fe do (20) 


em que o factor fora do integral depende da energia 
incidente E; e difundida Er, das constantes do electrão 
(cem) e das características do oscilador (massa redu- 
zida M, energia ho e momento du/ds). O integral é 
uma grandeza adimensional, função de E(;/E; e do 
ângulo de incidência O, sendo 


Hor 


f= Ea 
r 


K; 
Ext r p* 
2 J K, (r) 


“e (r) dir (21) 
a amplitude de difusão e O o ângulo sólido em torno 
da direcção especular. 

No cálculo das secções eficazes é possível separar 
as amplitudes de difusão f pef , associadas às projec- 
ções do momento dipolar respectivamente segundo uma 
direcção perpendicular e paralela à superfície [33]: 
E ————— K 

2 puta 2 
E + (Ki, Kr,) 
— 2K:, — Kr,) sen (Ki, -— Kr JZo | nes 


[4 


= ap eae 


E | q; eq Zo 
q +K HKD MO 
— 2(Ki iu Ke) sen (Ki; ai Ke pZo| (22) 


Ki 


Lo =D|/——---——————— -n £o 4 
A + (1, = | eq, ni 


a 
dy 
+ 2 cos(Ki;, — Ke )Zo | — 


Ki 4 | o 4 
= et À ed 
Ly E (&, + 1 
+ 2 cos(Ki, + Kr )Zo | (23) 


A secção eficaz diferencial 


do [EN1/2m 1º0,/dey 
dQ dE = E) M Eha Ea) Vero 


é proporcional a f?, que por sua vez é dado por 


E = IE, |? cos? O, + |f,|º sen? G, cos? £ (25) 


2 
M | 
em que & é o ângulo entre as componentes paralelas 
do momento dipolar, /*e do vector de onda q y defi- 
nido por 


l 
cos 8=— Ki, cos yo— Kr, cos (P'— qo)| (26) 


dy ” 


No caso de uma molécula adsorvida na superfície 
metálica, os electrões de condução do metal tendem 
a blindar o campo eléctrico do dipolo oscilante. Tra- 
ta-se de um efeito quase instantâneo quando a frequên- 
cia dos osciladores moleculares for muito inferior à 
frequência do plasmão dos electrões metálicos [36]. 
Este efeito de blindagem é usualmente descrito, num 
modelo clássico, com a introdução do conceito de 
dipolo imagem fictício, situado a uma distância, —Zo, 
abaixo da superfície [36, 37], cujas componentes, em 
conformidade com a fig. 7, são: 


dito, 
im — prim 
E 


Em termos deste sistema dipolo-dipolo imagem, as 
amplitudes de difusão totais serão: 


Re = É, (Zo) + f (—Zo) 
à —— É, (Zo) + t, (—Zo) 
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Fig. 8: Representação esquemática do dipolo |, e seu dipolo 
imagem E 
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Fig. 9: Variação da secção eficaz total em função do produto 
KZo para um electrão incidindo segundo O = 45º e 
com uma energia E, = 5,26 hw (a) Dipolo normal à 
superfície (b) Dipolo paralelo à superfície. Valores de 
E correspondem a osciladores de energia hwy = 0,25 eV 


situado a Zo = 2,5 A. Reprodução de [33 e 36]. 


A fig. 9 permite verificar o efeito da contribuição 
do dipolo imagem em termos do produto K Zo. Assim, 
quanto mais pequeno for este, maior será essa contri- 
buição. No caso limite de Zo = O a componente nor- 
mal de [' é 


E => f (Zo) 
enquanto que a componente paralela é nula. O dipolo 


imagem gera, nestas circunstâncias, a eliminação da 
interacção com os modos que oscilam paralelamente 


à superfície e aumenta de quatro vezes a secção eficaz 


de interacção com os modos perpendiculares. Isto sig: 
nifica que só estes devem ser observáveis, o que cor 
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responde à regra de selecção dipolar conhecida da 
espectroscopia do infravermelho. 

Para o Zo não nulo, mas muito inferior ao com- 
primento de onda associado ao electrão, a contribuição 
das vibrações normais à superfície mantém-se prepon- 
derante. Pelo contrário, quando À for da ordem de 
grandeza de Zo, a influência do dipolo imagem é for- 
temente atenuada e a interacção com os modos para- 
lelos devém importante (fig. 9). Isto constitui uma 
diferença considerável entre a espectroscopia vibracio- 
nal de perda de energia de electrões e a espectroscopia 
d> infravermelho. 

A modulação de É, e É, 
ções 22 e 23, permite prever que, modificando a ener- 
gia incidente, é possível variar a contribuição relativa 
dos modos perpendiculares e paralelos. Na fig. 9 
é visível que uma energia E; = 11,8 eV privilegia É, 
enquanto que t, é preponderante para E; = 3,65 . 
Nestas cdvcrenalicias. o estudo da evolução dos espec- 
tros em função da energia primária deve fornecer 
informações estruturais sobre a configuração das espé- 
cies adsorvidas [33]. 


explicitada nas equa- 


180º 0º 180º q 


B=45" 
ne 


hoo 


e — 
E e 
= mm cm 


50 


e = 
mm e a mm 


Precios 


ta 


E 
m 
= 
um 
- 


LES 
Re om 
ee É 


La 
mm 

Es mm 
q mm 


>» 
| 


[e 
“ 
q 
[o 
q 

É 


tc) 


Secção eficaz diferencial no caso de um dipolo na 
superfície (Zo = 0) para diferentes ângulos de inci- 
dência. Os resultados estão indicados para 4 =0 e 
9º = 0. Traço contínuo E, = 1,32 eV, tracejado E, — 
= 2,5 eV (Reprodução de [36]). 


Fig. 10: 


A fig. 10 mostra a influência dos ângulos de inci- 
dência, de análise e de orientação do dipolo no caso 
particular de Zo nulo. Nela são considerados três 
ângulos de incidência, próximo da normal (5º), quase 
rasante (85º) e intermédio (45º) para os casos de elec- 
trões primários de energia iguais a 2,5 e 1,32 eV. 
Exceptuando a incidência quase normal, em que a 
probabilidade de transmissão do feixe é elevada, a 
secção eficaz apresenta um máximo na vizinhança da 
direcção especular. Isto significa que, pequenos desvios 
no alinhamento do analisador em relação à direcção 
especular, podem provocar erros de medida, sobretudo 
se os resultados forem referenciados ao pico elástico. 
No entanto, tais desvios podem ser interessantes por 
atenuarem o fundo contínuo elástico e isso permitir 
uma melhor separação das perdas devidas às vibrações. 

Para Zo diferente de zero, a contribuição dos dipo- 
lo» é modulada pelo ângulo /&. Se a adsorção das molé- 
culas for aleatória, o termo cos'8 tem um valor médio 
1/2. Acontece, porém, ser frequente as moléculas 
adsorvidas tomarem orientações bem definidas em 
relação aos eixos privilegiados das superfícies. Nestas 
circunstâncias haverá valores particulares do ângulo 8, 
o que abre perspectivas quanto à possibilidade de se 
obterem, a partir destes espectros vibracionais, infor- 
mações sobre a orientação das moléculas adsorvidas 
em relação a uma superfície e mais concretamente em 
relação a eixos privilegiados desta última. 


2.2.2 — Interacção ressonante 


No estudo da interacção dipolar despreza-se o termo 
do potencial de curto alcance, admitindo que o tempo 
de permannécia dos electrões na zona de influência 
deste potencial é suficientemente pequeno para poder 
provocar uma excitação. Esta hipótese deixa de ser 
vélida sempre que se observar um mecanismo que 
obrigue o electrão a estar um tempo prolongado junto 
da molécula, Um desses processos consiste na captação 
do electrão num nível antiligante da espécie adsorvida, 
originando-se deste modo um ião negativo temporário. 
Na verdade, este mecanismo domina as secções efica- 
zes diferenciais vibracionais das moléculas em fase 
gasosa [11] e a sua secção eficaz aumenta significati- 
vamente quando a energia corresponder à energia de 
ressonância [38]. 

Um mecanismo de difusão por ressonância foi con- 
siderado por Ho [39] e uma primeira abordagem foi 
realizada por Davenport [40]. Segundo o tratamento 
deste autor, as propriedades da interacção ressonante 
têm que ver com as características do estado ligado 
nc qual o electrão fica temporariamente retido. 

Em fase adsorvida, em que a posição da molécula 
está fixa, a probabilidade de formar temporariamente 
um ião negativo não depende só da energia do elec- 
trão incidente, mas também da sua direcção, que deve 
ser adaptada à simetria do ião formado. A distribuição 


angular dos electrões reemitidos reflecte, por outro 
lado, a simetria do estado ressonante, Estas duas pro- 
priedades podem pois ser utilizadas para determinar a 
geometria das ligações da molécula com o subtsracto. 

O tratamento de Davenport mostra que a difusão 
ressonante não implica nenhuma regra de selecção 
(permitindo observar vibrações paralelas à superfície) 
e não produz unicamente nas direcções vizinhas da 
especular. Do ponto de vista experimental, isto deveria 
permitir distinguir a difusão ressonante da interacção. 
Acontece, porém, que o tempo de vida dos iões nega- 
tivos é muito pequeno em fase adsorvida, o que difi- 
culta a observação da difusão ressonante, 


3 — DISPOSITIVOS EXPERIMENTAIS 


Os aspectos teóricos relacionados com a espectro- 
metria de perda de energia de electrões que referimos, 
definem, como grandezas físicas importantes nesta téc- 
nica, a energia E;, dos electrões incidentes, a ener- 
gia, 4 E, cedida por estes às moléculas do meio difusor 
e os ângulos de incidência, O, e de difusão 9". Agindo 
nestas grandezas é possível programar três formas de 
actuar experimentalmente. 


— a energia e a direcção do feixe incidente, bem 
como a direcção de difusão são mantidas cons- 
tantes e analisa-se em energia o feixe difundido, 
pondo-se em evidência as diferentes perdas de 
energia, 

— as direcções dos feixes e um valor de perda de 
energia são fixos variando-se a energia incidente 
(recolhem-se assim, informações sobre o meca- 
nismo de excitacçção da perda 4 E escolhida). 

— Mantém-se constante a energia, o ângulo de inci- 
dência e a perda AE e faz-se uma análise da 
distribuição angular dos electrões difundidos. 


Para permitir realizar estes tipos de experiências, 
um espectrómetro terá que contar as seguintes três 
unidades, mais ou menos sofisticadas: 


— Um monocromador englobando uma fonte de 
electrões (cátodo mais canhão a electrões), um 
selector de velocidades e uma Óptica de saída. 

— Uma câmara de difusão onde se coloca a amos- 
tra ou se injecta um gás. 

— Um analisador de electrões compreendendo uma 
óptica de entrada, um selector de velocidade e 
um detector de electrões. 


Num dos primeiros espectrómetros, montado por 
Ramsauer [41] o conjunto destas unidades era colo- 
cado num campo magnético uniforme. Tratava-se de 
um aparelho simples, cujo funcionamento se baseava 
na não detecção, por uma gaiola de Faraday, 
de todos os electrões que sofriam uma colisão com 
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Fig. 11: Esgema simplificado do dispositivo de Schulz (Repro- 
dução de [43]). 


transferência de momento não nula. Permitia estudar a 
variação da intensidade de corrente eléctrica na gaiola 
de Faraday, em função da intensidade do campo mag- 
nético e por isso só fornecia dados sobre secções efi- 
cazes globais. 

Shulz [42] desenvolveu um aparelho cujo princí- 
pio ce funcionamento consistia em captar os electrões 
que, após uma colisão, ficavam com uma energia 
E'< e Vp, sendo V, a profundidade de um fosso de 
pote: cial criado na câmara de colisão. Um esquema 
sucinto deste espectrómetro está representado na 
fig. 211 [43]. 

O feixe de partículas, emitido por um canhão de 
electrões, é focalizado por um campo magnéstico coa- 
xial, de fraca intensidade (100 a 500 Gauss) e acele- 
rado por uma diferença de potencial, Va, aplicada 
entre F e a câmara de colisão. Esta, é um cilindro, 
centrado no eixo do feixe e constituída por uma gre- 
lha G, limitada, nas suas extremidades, por duas placas 
paralelas dotadas de um diafragma circular. Uma dife- 
rença de potencial aplicada entre M e G, cria um fosso 
de potencial Vp (= 1% da tensão entre M e G) ao 
longo do eixo do cilindro. Todos os electrões de energia 
residual E” < e V, ficam retidos no eixo e difundem 
posteriormente até ao colector M. As restantes parti- 
culas são recolhidas em N. 

Medindo a variação de intensidade de corrente no 
circuito do eléctrodo M, em função da energia dos 
electrões incidentes, observam-se máximos da curva 
sempre que a energia E” estiver compreendida entre 
E, e E, + eVp, sendo E, a energia de um nível exci- 
tado. 

O espectrómetro de Shulz adapta-se bem ao estudo 
de transições com grande transferência de momento 
não descritos pela aproximação de Born, pois o seu 
funcionamento parte da detecção de electrões que cede- 
ram ao alvo praticamente toda a sua energia cinética. 
É, por outro lado, dotado de uma grande sensibilidade 
pois, a coleta de electrões pelo eléctrodo M, tem um 
rendimento muito perto dos 100 %. Um inconveniente 
deste aparelho é o de não fornecer dados sobre a 
difusão angular, uma vez que a detecção de electrões 


= 


é feita segundo uma geometria 4 7. 


Este tipo de dispositivo foi retomado por outros 
autores [11, 44, 45], mas mais recentemente Sanche 
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[46, 47] adaptou-o de forma a permitir a análise, em 
transmissão, de substâncias sólidas. Os electrões orien- 
tados por um campo magnético são seleccionados por 
um monocromador trocoídeo e incidem na amostra. 
Esta é depositada, sobre uma fita de tungsténio arre- 
fecida à temperatura do azoto líquido, por condensação 
de um vapor. A intensidade de corrente devida aos 
electrões que atravessam o alvo é então medida em 
função da energia incidente. 


A selecção dos electrões difundidos pode ser feita 
com o auxílio de um potencial retardador. É o que se 
passa no espectrómetro sugerido por Kupperman [48, 
49]. Uma representação esquemática permitindo com- 
preender o princípio do seu funcionamento está indi- 
cada na fig. 12. 

O canhão de electrões compreende um cátodo e 
um diagrama F levado a uma tensão ligeiramente nega- 
tiva em relação ao cátodo. 

A câmara de colisão é limitada por um eléctrodo 
plano com um diafragma circular e por uma grelha G,, 
cilíndrica. O conjunto é polarizado positivamente em 
relação a F, de forma a acelerar os electrões incidentes. 


O analisador de electrões difundidos é constituído 
por um conjunto de grelhas, G> e Gs, e por um eléc- 
trodo colector, M. Os electrões não difundidos são 


detectados por uma gaiola de Faraday, diante da qual 
uma grelha, G”, se destina a repelir os electrões secun- 


dários arrancados às paredes da gaiola. 

Para uma dada energia dos electrões incidentes, 
polariza-se a grelha Gs, de forma a que nenhum elec- 
trão difundido possa atingir o eléctrodo M. Diminuindo 
progressivamente o campo retardador existente entre G; 
e G>, vamos detectando os electrões difundidos come- 
cando portanto, pelos mais energéticos. Recolhem-se, 
deste modo, curvas integrais de intensidade de corrente 
em função do potencial retardador. As curvas terão 
então que ser diferenciadas para se obterem os espec- 
tros. 

Este tipo de dispositivo não fornece também infor- 
mações angulares e a existência de uma corrente resi- 


| 
| me ita 


E E —— = Ts E — 


Fig. 12: Esquema simplificado do espectrómetro de Kupper- 
mann (Reprodução de [43]). 


dual limita o valor mínimo de energia com que se 
podem efectuar medidas. 

O esquema da fig. 12 refere-se a experiências com 
alvos gasosos, mas por uma mudança de geometria é 
possível trabalhar com amostras condensadas [50]. 

O estudo experimental da influência dos ângulos 
de incidência e de difusão nas secções eficazes de exci- 
tação, só é possível se o monocromador e o analisador 
e electrões puderem rodar um em relação ao outro. 
O primeiro dispositivo com esta característica foi apre- 
sentado no início dos anos 30 por Bullard [51]. Pos- 
teriormente, assistiu-se a um desenvolvimento enorme 
deste tipo de espectrómetros de electrões. Klein [43] 
historia essa evolução até 1968 e dos trabalhos por ele 


referidos citamos, pela sua importância, os de Hugher 
[52] e Ehrhardt [53, 54]. Depois de 1968, são de 
referir as contribuições de Ballu [55, 56], Read [56, 
58], Preston [59], Roy [60], Froitzheim [61], e 
Ergenton [62, 63.] Os artigos de revisão de Celotta 
[15] e de Raether [64] contêm citações com interesse. 

Neste artigo desenvolvemos, com algum detalhe, 
um dispositivo, semelhante ao de Roy, montado por 
Rei Vilar [7, 65], e especialmente destinado à obten- 
ção de espectros de perda de energia de electrões len- 
tos rectrodifundidos por alvos condensados. A repre- 
sentação esquemática deste espectrómetro está indicada 
na fig. 13. 

Os electrões emitidos por um filamento de tungs- 
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Fig. 13: Esquema de um espectrómetro de electrões lentos com possibilidade de análise angular [7, 60]. 
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ténio revestido a tório, são acelerados a cerca de 6 eV 
e entram através de um diafragma num monocromador 
electrostático, C;, de forma cilíndrica (27%). O campo 
eléctrico é mantido entre duas grelhas, G, e G>, e os 
electrões não seleccionados são capturados pelas pla- 
cas de cobre P, e P>. A energia dos electrões incidentes 
pode variar entre o e algumas centenas de eV e é 
definida pelo potencial aplicado à amostra. Esta é 
colocada na câmara de ionização mantida equipoten- 
cial. Os electrões difundidos são seleccionados no ana- 
lisador, G>, idêntico ao monocromador e detectados 
por um multiplicador de electrões (Chaneltron). Os 
impulsos são amplificados, discriminados, contados e 
recebidos num sistema de análise de sinais (analisador 
multinacional e/ou microcomputador). Um sistema de 
variação contínua de tensões (rampa) aplicado, por 
exemplo, ao analisador, varre a gama de perda de 
energia dos electrões, desde zero (pico elástico) até 
à energia incidente (secundários verdadeiros) o número 
do canal do analisador multicanal. Todo o sistema é 
blindado electrostáctica e magneticamente e fica sob 
um vácuo de = 10º Torr. 

Terminamos este resumo sobre dispositivos expe- 
rimentais, assinalando a associação que muitos auto- 
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Fig. 14: Espectro de excitação do benzeno (Reprodução 
de [68]). 


res têm feito de um espectrómetro de electrões a um 
microscópico electrónico. Um trabalho de revisão de 
progressos obtidos neste campo foi publicado por Col- 
liex [66]. No entanto para uma visão do muito que 
se tem publicado sobre este assunto, convidamos o 
leitor interessado a consultar um artigo recente de 
Egerton [67], onde este autor compilou cerca de 
700 referências. 
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Fig. 15: Espectro de perda de energia de electrões em benzeno 
obtido para várias energias incidentes (a) 6 eV (b) 
&eV (c) lleVe (d) 14 eV (Reprodução [69]). 
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4 — EXEMPLOS DE ESPECTROS DE PERDA DE 
ENERGIA DE ELECTRÕES 


Reproduzimos neste parágrafo alguns espectros de 
electrões difundidos, procurando dar uma ideia das 
potencialidades da técnica que temos vindo a analisar. 
A selecção que fizemos é necessariamente limitada 
pois, são muitos os resultados experimentais até agora 
publicados. Assim, concentrámos a nossa atenção 
essencialmente sobre três questões: a excitação de 
níveis de energia dificilmente acessíveis por meios ópti- 
cos, a possibilidade de recolher certas características 
moleculares através de análise adequada dos espectros 
e o estudo de substância adsorvidas em superfícies 
metálicas. No final do parágrafo indicaremos algumas 
referências sobre outros tipos de aplicações eventuais 
desta espectroscopia. 


4.1 — Estudo das excitações moleculares do benzeno 


Para mostrar o interesse da espectroscopia de perda 
de energia de electrões na identificação de estados de 
energia, não observáveis facilment por espectroscopia 
óptica, optámos por seleccionar trabalhos publicados 
sobre uma substância: benzeno. 

O espectro da fig. 14, publicado por Compton [68] 
e obtido por transmissão através de filmes finos de 
benzeno, dá uma ideia do tipo de informação reco- 
lhida. Os picos obtidos identificam-se bem com os 
níveis energéticos tripleto e singuleto do benzeno, o 
que mostra a possibilidade da técnica no estudo das 
transições opticamente proibidas. O pico intenso de 
perda de energia, observado a 1,55 eV, está abaixo 
do estado electrónico e foi interpretado como corres- 
pondendo a um ião negativo de ressonância. 

Utilizando um espectrómetro do tipo Kuppermam, 
Merkel [69], estudou a evolução dos espectros em 
função da energia dos electrões incidentes (6, 8, 11 e 
14 eV) fig 15. É possível verificar que a excitação dos 
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Fig 16: Secção eficaz de excitação dos níveis T,, S, e S, do 
benzeno, em função da energia incidentes (Reprodu- 
ção de [69]). 
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Fig. 17: Espectro de perda de energia de electrões de 13,6 eV 
em benzeno, para valores de 0'=9º, 23º e 54º 
(Reprodução de [70]). 
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níveis tripletos é mais importante para energias pri- 
márias mais baixas (= 10 eV). A fig. 16 mostra valo- 
res relativos das secções eficazes de excitação para os 
estados T; (interdita por spin), S; (interdita por razões 
de simetria) e S: (opticamente permitida). 

A influência da direcção de difusão sobre os espec- 
tros foi estudada por Doering a partir de espectros 
obtidos com uma energia incidente fixa (13,6 eV) mas 
detectando os electrões difundidos, segundo três ângu- 
los diferentes (0º = 9º, 22º e 54º). Os resultados deste 
autor (fig. 17) confirmam a importância do ângulo 6º 
na observação das excitações tripletos. Na verdade, a 
probabilidade de registo de uma transição singuleto-tri- 
pleto aumenta para valores crescentes do ângulo de 
difusão, como aliás, se depreende melhor da fig. 18. 

Jogando simultaneamente com a energia incidente e 
o ângulo 6”, podemos portanto encontrar as melhores 
condições experimentais de análise de uma dada exci- 
tação. A fig. 19 representa a zona do espectro do 
benzeno correspondente ao intervalo de perdas com- 
preendido entre 4,3 e 5,5 eV (inclui excitações T; e 
S1), obtida em três condições experimentais diferentes 
(três pares de valores E;(09, [71]. Observa-se deste 


21 


Benzeno 

E. = 13,6 
1 

0,007 Torr 


—0-9O ; DO) em (E | 


10º 


10" 
ad ae am 

| q E 
E: 28 Ri 
E , e 
= 10 Ao 
E 8 
a é O 5/53 
E a 
z nd OIT; +9/)/53 

“ 

o Ts 


ie 


o O ao 60º Bo 


Fig. 18: Razões de intensidades de várias transições observa- 
das no espectro do benzeno em função do ângulo de 
difusão 6'. Energia incidente 13,6 eV (Reprodução 
de [70]). 


modo que E; elevados e O” pequenos (parte a da fig.), 
permitem individualizar a transição So” S,, incluindo 
excitações vibracionais do estado S;, do benzeno. Pelo 
contrário, valores baixos de E; e elevados de O”, favo- 
recem a observação da transição interdita So” T; 
(parte c da figura). Nas condições intermédias da 
fiz 19b, o espectro é menos explícito devido à mis- 
tura das duas excitações. 

Christophorou, no seu artigo de revisão [27], 
ocupou-se do pico que aparece a 1,55 eV na fig. 14, 
relacionando-o com a formação de um ião negativo 
temporário. São assim, referidos trabalhos de Boness 
[72] e de Sanche [73] mostrando a estrutura vibra- 
cional do pico correspondente ao ião de ressonância 
(ver fig. 20) e de Larkin [74], que interpretou a 
natureza do processo envolvido. 

Outros resultados experimentais sobre o benzeno 
poderiam ainda ser citados. Pensamos porém, que os 
indicados são susceptíveis de mostrar a importância da 
espectroscopia de perda de energia de electrões na 
observação e identificação da natureza de excitações 
moleculares. 


4.2 — Espectros de perda de energia e teoria dieléc- 

trica 

Fazemos agora uma breve referência à interligação 
entre os espectros de perda de energia e a teoria die- 
léctrica. Conforme referimos no ponto 2, é possível 
reproduzir um espectro a partir do conhecimento da 
função de perda I(— 1/£) calculada a partir do conhe- 
cimento da constante dieléctrica do meio (ver fig. 4). 
Assinalámos igualmente ser possível o caminho inverso: 
a partir do espectro determinar I(— 1/€), a constante 
dieléctrica e finalmente certas características molecula- 
res. Daniels [22] e Raether [64] dão vários exemplos 
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de trabalhos deste tipo. Neste trabalho, fazemos refe- 
rência a contribuições de Ritsko para o estudo do 
poliestireno e do poliacetileno. 

No caso do poliestireno [25], a função de perda 
fo' calculada a partir de um espectro de excitação de 
valência e do cor de electrões, estendendo-se num 
intervalo de perda de energia de 1 a 400 eV e irra- 
diando com electrões de 80 KeV. A partir da função 
de resposta foi caclculada a reflectividade e o índice 
de refracção, tendo-se verificado um acordo excelente 
com resultados experimentais obtidos por via óptica. 
É o que se pode observar na fig. 21 no que respeita à 
reflectividade. Nesta figura indicam-se ainda coeficien- 
tes de absorção óptica calculados a partir da função 
de perda. 

Para estudar a estrutura electrónica do poliaceti- 
leno Ritsko [75, 76] analisou os espectros de filmes 
finos desta substância, pura e contendo AsF; ou 1; 
como impurezas. 
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Fig. 19: Região de perdas A E situada entre 4,3 e 5,5 eV de 
espectros de perda de energia de electrões em ben- 
zeno, obtidos para três pares de valores Eq”: 
(3ileV — 10º), (184eV-— 16%) ce (12533€eV — 30. 
(Reprodução de [27]). 


INTENSIDADE (U.A.) 


INTENS IDADE 


os LO 1.2 1,4 1.6 1.8 
ENERGIA E, (eV) 


BF lmevi: faiz ef 110 | 112 -|4 1064 
FG I 2 3 a 


Ss 4 Lisias. ros igGÃãl 
hE (ev) 


Fig. 22: Espectros de perda de energia de electrões de 80 KeV 
em poliacetileno puro (.), dopado com 1, ([]) e com 
AsF; (A). Reprodução de [76]). 
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Como se verifica na fig. 22, o espectro do poliace- 
tileno puro é caracterizado por um limiar de absor- 
ção de 1,4 eV e um pico alargado centrado em 4,5 eV. 
As curvas referentes ao polímero dopado evidenciam 
um aumento de absorção na zona do «gap», assim 
como um novo pico à volta de 1 eV, de fraca reso- 
lução, mas mais pronuciado no caso do iodo. A figura 
mostra ainda, que a introdução deste último elemento 
provocou uma mistura mais heterogénea de regiões 
dopadas e não dopadas do que a do AsF; (o pico a 
45 eV persiste ainda). 

Aplicando a análise de Kramers-Kronig [22], Ritsko 
calculou as curvas referente às partes real e imaginária 
da constante dieléctrica indicadas na fig. 23. Reco- 
lhendo os valores de e, e de e; da parte a da figura e 
calculando Es (e + e), é possível observar um máximo 
a 1,8 eV na função de perda que foi interpretado como 
correspondendo a uma absorção interbanda. A parte b 
da figura evidencia que o pico intenso de e; se deslo- 
cou, por efeito da dopagem, para as energias muito 
baixas, o que é consistente com o desaparecimento do 
gap característico de um comportamento semicondutor. 
As excitações dominantes passam a ser próprias de um 
metal. A análise de Kramers-Kronig mostra que o pico 
oservado a 4,1 eV ocorre depois de £, passar por zero, 
numa região onde £, é relativamente pequena e decresce 
rapidamente (fig. 23a). Este pico é então interpretado 
como um plasmão associado com excitações colectivas 
dos electrões 7. 


4.3 — Análise da emissão secundária 


A análise da parte do espectro correspondente à 
emissão de electrões secundários, apresenta, para algu- 
mas substância, ligeiras estruturas, igualmente observa- 
das em experiências de fotoelectrão [77-80]. Como se 
pode verificar pela figura 24, a posição em energia das 
estruturas é independente da energia da partícula inci- 
dente e da natureza da radiação (electrões ou fotões), 
o que parece indicar que a sua origem é independente 
destes factores. Mecanismos susceptíveis de contribui- 
rem para o aparecimento destas estruturas podem em 
princípio ligar-se a: 


— excitação de electrões para a banda de condução 

— difusão inelástica de electrões incidentes 

— electrões excitados por declínio Auger ou autoio- 
nização de excitões de energia elevada 


A interpretação de Ueno dos resultados da fig. 24 
[79] não contempla estes Últimos mecanismos. O autor 
mostra como, para as substâncias consideradas, alca- 
nos, a análise dos «acidentes» da emissão secundária 
permite tirar conclusões sobre a banda de condução, 
nomeadamente qua sua estrutura é quase independente 
do número de átomos de carbono. 

A autoionização de excitações de energia elevada é 
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Fig. 24: Estruturas da emissão de electrões secundários: (1) 
Espectro de perda de energia de electrões de 20 eV 
em: (a) n-C,aHsg, (b) n-CyHo,, (Cc) n-CyHoo, (d) polie- 
tileno; Espectro de fotoelectrão U.V. (21,2 eV) em: 
(e) n-CHsg, (É) n-CyHhy, (h) — polietileno: II Espec- 
tro de perda de energia de electrões em nCyH.. 
(Reprodução de [79]). 
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Fig. 25: Espectro de electrões secundários de perileno policris- 
talino para várias energias incidentes. Abcissa infe- 
rior: perda de energia, Abcissa superior: energia 
cinética E, dos electrões secundários. (Reprodução 
de [77]). 
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Fig. 26: Espectro de perda de energia de electrões de 34 eV 
em poliestireno. 


observável noutros trabalhos de Ueno [77, 78] sobre 
cristais orgânicos. É o caso das estruturas À e B atri- 
buídas aos electrões emitidos no declíneo das excita- 
ções correspondentes às perdas de energia D (6,42 eV) 
e E (8,3 eV). Tomando como energia de ionização do 
pcrileno o valor de 5,8 eV obter-se-ia, para a energia 
cinética dos electrões correspondentes às estruturas A e 
B, respectivamente 0,62 e 2,5 eV, valores não muito 
diferentes dos determinados por espectroscopia do 
fotoelectrão [81]. A natureza da bossa centrada à volta 
da perda de energia de = 5 eV foi interpretada como 
consequência da estrutura da banda de condução [78]. 


Um resultado semelhante ao obtido por Ueno com 
cristais, foi observado na emissão secundária em fil- 
mes de poliestireno [7]. Com efeito, parece existir uma 
correlação entre a banda centrada a 13 eV e identifi- 
cada como um plasmão característico do poliestireno 
[25], e o cotovelo que se apercebe à volta de 29 eV 
no pico dos secundários (fig. 26). Na verdade, admi- 
tindo que a energia de um estado superexcitado é 
igual à do plasmão que lhe deu origem, podemos cal- 
cular, a partir da figura, para a energia de ionização 
do poliestireno um valor = 8 eV, que está de acordo 
com dados da literatura [82, 83]. 


4.4 — Espectros vibracionais de substâncias adsorvidas 
em superfícies metálicas 


Os exemplos que vamos citar foram seleccionados 
dos vários apresentados por Bertolini [33] e Lopez 
[34]. A figura 27 refere-se aos espectros obtidos por 
Barnes [84] para estudar a adsorção de hidrogénio 
sobre a face (100) do tungsténio. Nas partes a e c da 
da figura, o feixe de electrões difundidos é observado 
segundo a direcção especular, 4 O = O, enquanto que 
bed foram medidas fora dessa direcção, 40 = 25º, 
Para uma dada direcção de difusão os espectros dife- 
rem entre si, por um deles corresponder ao caso da 
superfície estar saturada com hidrogénio. Em todas as 
experiências a energia incidente foi de 55 eV e o 
ângulo de incidência 60º. 

No caso da superfície saturada em hidrogénio 
(Ou = 2) a análise dos electrões rectrodifundidos na 
fig. 27a revela um só pico a 130 mev, que deve 
ser associado a um modo de vibração, v;, perpen- 
dicular à superfície (recorde-se que os modos nor- 
mais têm uma contribuição preponderante). Por 
outro lado, a existência de um só modo leva a 
supor que o hidrogénio adsorvido deve estar em posi- 
ção linear (átomos de H por cima de átomos de W) 
segundo uma simetria C ,) de ligação ao substrato. 

A recolha do espectro fora da direcção especular 
faz aparecer duas novas vibrações (» a 80 meV ev; a 
160 meV), que são identificadas com vibrações de um 
hidrogénio posicionado como se indica na fig. 27b. 
A geometria das ligações do átomo de hidrogénio é 
então em ponte, com simetria C»v. 

Comparando os espectros a e b com os indicados 
em ce d, observa-se que as energias de vibração variam 
em função do teor de recobrimento em hidrogénio. 
Nota-se uma passagem de uma estrutura ordenada 
C(2 x 2)H se Oy < 0,5, para a estrutura p(1 x DH 
se Oy — 2 e que ela está ligada às variações de energia: 


C(2x2)H 
ptix DH 


v,i— 155 meV 
V|— 150 meV 


v;— 6O meV 
v,— 80 meV 


vY— 125 meV 
vw— 160 meV 


Continuando a observar 3 vibrações, isso significa 
que a ligação hidrogénio continua a ser em ponte. 
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A modificação do valor das energias de vibrações tem 
que ver com o valor do ângulo « definido por W-H-W 
que é inferior a r/2 no caso da estrutura C(2 x 2) e 
superior a este valor para a estrutura p(l X 1). 

Um outro exemplo com interesse, corresponde à 
adsorção do CO. A fig. 28 reproduz os picos caracte- 
rísticos dos espectros vibracionais de CO adsorvido 
em Pt (1,1,1) [85] e Ni(1,1,1) [33]. Nas experiên- 
cias sobre platina e com um recobrimento O < 0,33 só 
se observa uma espécie de CO adsorvida, caracterizada 
pela vibração de valência C-O(vco— 259 mev) e C-me- 
tal (veu — 58 meV), fig. 28a. Para um recobrimento 
tal que 0,33 < Oco < 0,5 (estrutura C(4 X 2)), nota-se 
a presença de duas espécies de CO adsorvidas, fig. 28b, 
caracterizadas por vco — 126 meV e vem — 58 meV 
(semelhante à anterior) e por vco — 232 meV e vem — 
— 45 meV. No caso de adsorção sobre Ni(ll1l1) a 
forma do espectro mantém-se quando variam as razões 
de recobrimento em CO, fig. 28c, e assemelha-se à 
estrutura que se obtém para 0,33 < Oco < 0,5 no caso 
da Pr (111). Assinala-se a presença de um pico de 
fraca intensidade a 79 meV, associado a uma vibra- 
ção vem assimétrica (paralela à superfície). 


A observação de um só modo vibracional C-metal 
(58 meV) na adsorção sobre Pr com baixo teor de 
recobrimento, sugere a existência de espécies lineares. 
Aumentando o recobrimento, é de esperar o apareci- 
mento de ligações em ponte. No níquel só é de esperar 
uma ligação em ponte. A energia da vibração Vco é 


INTENSIDADE 


muito sensível ao local de adsorção. O valor 260 meV 
está associado ao CO linear enquanto que os 230 meV 
ac CO em ponte. O estudo da Quimisorção do CO 
pode deste modo fornecer dados sobre os locais de 
adsorção. 

Como último exemplo, vamos referir a adsorção 
de hidrocarbonetos em C; sobre platina e níquel. 
À observação dos espectros de vibração de hidrocar- 
bonetos em C;>, adsorvidos sobre platina, em diferentes 
condições experimentais, fig. 29, mostra a sensibilidade 
da técnica na observação da natureza dos grupos que 
aparecem à superfície. 


A parte a da figura 29 refere-se a experiências a 
uma temperatura abaixo de 300 K, com C;H; e C.D:» 
adsorvido. O pico a 162 meV, não é afectado pela 
deuteração, pelo que deve corresponder ao modo de 
elongação vc c. Os outros picos implicam movimentos 
de átomos de hidrogénio e foram atribuídos a vibra- 
ções de deformação (95 e 122 meV) e de elongação 
(373 meV) de C-H. A grande intensidade do pico a 
95 meV leva a pensar que ele se refere a um modo 
normal à suprfície, estando os outros principalmente 
polarizados paralelamente à superfície. 

A temperaturas inferiores a 260 K, a adsorção de 
etileno, com teor de recobrimento elevado, conduz ao 
espectro da fig. 29b. Os picos a 176, 123 e 98 meV 
identificam-se com grupos CH». A vibração a 364 meV 
revelaria uma hibridação do tipo sp' para os átomos 
de carbono. A perda a 152 meV corresponde a um 
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Fig. 27: Espectros vibracionais de hidrogénio adsorvido em 
W(100) obtidos com electrões de 5,5 eV incidindo 
segundo 6 = 60º (Reprodução de [33]). 
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modo de vibração ve-c polarizado paralelamente à 
superfície. 


Os espectros da parte c da fig. (29) foram obtidos 
a uma temperatura superior à ambiente, com etileno 
adsorvido, ou após hidrogenação do acetileno adsor- 
vido. A comparação destes com os anteriores evidencia 
uma modificação profunda das espécies adsorvidas. 


O modo de vibração independente da deuteração 
situa-se agora a 140 meV. A intensidade elevada do 
pico desta vibração C-C leva a supor uma orientação 
normal à superfície da ligação C-C. Isto significa que 
se deve ter realizado uma reacção em superfície, acom- 
panhada da formação de uma espécie superficial está- 
vel. 

Como se depreende dos exemplos citados, é possí- 
vel observar directamente frequências de modos vibra- 
cionais em fase adsorvida, através da espectrometria 
de perda de energia de electrões lentos. Este facto, 
torna esta técnica útil na identificação da natureza dos 
complexos superficiais formados por quimidsorção, do 
número de ligação adsorvido-metal e da natureza dos 
locais de adsorção, etc. 

Os exemplos que citámos dão uma ideia das capa- 
cidades da espectroscopia de perda de energia de 
electrões. Como eles abrangem apenas uma pequena 
gama de substâncias é conveniente referir neste 
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Fig. 28: Picos característicos dos picos vibracionais de CO 
adsorvido em P, (111) e N, (111). (Reprodução 
de [33]). 


momento alguns estudos envolvendo outros tipos de 
meios. Assim, assinalamos trabalhos sobre gases raros 
[15, 22, 60], metais alcalinos [15] e de transição 
[22, 641], semicondutores [22, 88], sais de halogéneos 
[12, 22, 86, 87] e óxidos [15, 22, 89]. Citamos, 
ainda o artigo de Raether [64] onde o autor apresenta, 
em forma de tabela, uma lista de referências a traba- 
lhos envolvendo uma gama variada de elementos e 
compostos. Há, por fim, que realçar o facto da espec- 
troscopia por impacte de electrões poder aplicar-se em 
microanálise desde que seja possível dispor de uma 
biblioteca de espectros de perda de energia de subs- 
tâncias puras. É o que se depreende de trabalhos publi- 
cados sobre a análise de misturas gasosas [15] ou de 
amostras sólidas [90-92]. Os resultados obtidos com 
feixes de electrões, de energia suficientemente elevada 
para excitarem camadas internas, mostram que esta 
espectroscopia é particularmente interessante na detec- 
ção e análise de elementos leves. Nestes casos a sua 
sensibilidade é mesmo substancialmente maior do que 
a da microanálise por raios X [91]. A gama de per- 
das de energia explorável para fins quantitativos é, no 
entanto, limitada. Com efeito para valores de AE 
< 100 eV os picos correspondentes a excitações de 
camadas internas, podem ser submergidos num fundo 
motivado pela difusão inelástica nos electrões de valên- 
cia. Por outro lado, para perdas de energia muito ele- 
vadas (> 2000 eV), a intensidade dos electrões difun- 
didos é baixa misturando-se com o ruído da aparelha- 
gem. Este tipo de problemas tem sido abordado em 
vários artigos, em que se descrevem métodos de trata- 
mento dos espectros [65, 92-99] de forma a eliminar 
a contribuição do fundo e melhor evidenciar os picos 
que se pretende observar. 


5 — CONCLUSÃO 


Neste trabalho apresentámos alguns aspectos teó- 
ricos e experimentais característicos da espectroscopia 
de perda da energia de electrões difundidos inelastica- 
mnte. Trata-se de uma espectroscopia complementar 
de outras melhor conhecidas, distinguindo-se-lhes delas 
por, num mesmo dispositivo, ser viável explorar uma 
larga gama contínua de energias (1 a 10º eV) e obser- 
var perdas de energia extendendo-se desde as dezenas 
dc meV às centenas de eV. 

Os electrões incidentes podem classificar-se em 
três grupos consoante a sua energia cinética Es: Elec- 
trões lentos (< 100 eV) aos quais não é aplicável a 
aproximação de Born, Electrões rápidos (E, > 100 eV) 
passíveis de se lhe aplicar a aproximação óptica e pos- 
sibilitando a excitação de electrões do cor e Electrões 
com energia intermédia e por isso talvez, menos inte- 
ressantes. 

De entre os electrões lentos há que destacar os de 
energia inferior à dezena de eV, por permitirem obser- 
var perdas de energia da ordem das dezenas ou das 
centenas do meV, isto é. detectar e identificar modos 
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vibracionais de superfície. Nestas circunstâncias pode 
definir-se uma nova espectroscopia vibracional, par- 
ticularmente bem adaptada ao estudo configura- 
ciona! dos efeitos químicos estáveis, fixos à superfície 
de cristais metálicos. Torna-se por exemplo viável 
identificar complexos superficiais formados por qui- 
misorção, definir o número de ligações adsorvido- 
-metal e o local de adsorção, apreciar modificações de 
ligações químicas no interior das moléculas adsorvidas 
e consecutivas à quimisorção, etc. Um melhor aprovei- 
tamento da técnica neste domínio de perdas de energia, 
implica, no entanto, o desenvolvimento do estudo da 
interacção electrão-vibração. Uma outra particulari- 
dade dos electrões lentos, está no facto de induzirem 
a formação de plasmões de superfície, e excitações 
electrónicas de valência, não necessariamente concor- 
dantes com as regras de selecção óptica. Transições 
opticamente interditas por razões de spin ou de sime- 
tria são reveláveis nos espectros, dentro de uma gama 
de perdas de energia da dezena de eV. O decaímento 
dos plasmões de algumas dezenas de eV, pode ser 
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seguido através da observação de estruturas, na parte 
do espectro referente à emissão secundária. Enfim, 
analisando perdas de energia entre 1 e 50 eY, é possível 
recolher informações complementares sobre a estrutura 
das substâncias. 


Os electrões rápidos são especialmente interessan- 
tes, por provocarem a excitação dos electrões do cor 
e a observação de transições interbandas nos semicon- 
dutores. São características que permitem a utilização, 
da espectroscopia de difusão de electrões, em microaná- 
lise. Por outro lado, tendo em conta a aplicabilidade 
da teoria dieléctrica, é possível retirar dos espectros 
valores de constantes moleculares com interesse (absor- 
tividade, reflectividade, forças de oscilador ópticos). 


A gama de energias incidentes intermédias não 
parece tão interessante. As excitações por elas produ- 
zidas seguem as leis da selecção óptica sobretudo 
quando analisadas segundo pequenos ângulos de difu- 
são. À sua observação pode no entanto ser perturbada 
pela detecção dos electrões secundários. 


Fig 29: Espectros vibracionais de espécies estáveis formados 
sobre P, (111). (Reprodução de [33]). 
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RESUMO 


O projecto de circuitos com dispositivos semicon- 
dutores de potência requer uma caracterização física 
correcta do seu desempenho; para tal, torna-se neces- 
sário desenvolver modelos andliticos de fácil inserção 
em simuladores de circuito. Neste artigo, reunem-se 
as principais contribuições científicas do Prof. Pedro 
Teixeira na modelação e simulação de transistores 
bipolares e MOS de potência. Os exemplos apresenta- 
dos salientam a precisão dos modelos, a sua eficiência 
computacional e a sua utilidade para o projecto de 
circuitos. 


ABSTRACT 


The design of electronic circuits with semicon- 
ductor power devices quests for an accurate physical 
characterization of them performance. For this, analytic 
models that can be easily implemented in circuit simu- 
lators have to be derived. In this paper, the main 
scientific contributions of Prof. Pedro Teixeira in the 
area of bipolar and MOS transistors modelling and 
simulation are collected. Examples are presented that 
ascertain the accuracy of the models, their computa- 
tional efficiency and their usefulness in circuit design. 


NOTA PRÉVIA 


Com esta publicação, dá-se cumprimento ao desejo 
manifestado pela Secção de Electrónica e pelo Depar- 
tamento de Engenharia Electrotécnica e de Computa- 
dores do I. S. Técnico em assinalar de forma perdurá- 
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vel a importante obra científica do Prof. Pedro L. B. 
Teixeira. 


1. INTRODUÇÃO 


O projecto de circuitos electrónicos de potência 
com transistores bipolares de potência (TBP) torna 
necessária a simulação em computador do seu compor- 
tamento, para diferentes regimes de operação. A utili- 
zação de modelos numéricos para os TBPs exige usual- 
mente custos computacionais elevados. Além disso, 
requer o conhecimento completo, físico e dimensional, 
do dispositivo, o que por sua vez envolve um programa 
de medidas extenso para dispositivos existentes, ou 
uma simulação precisa dos processos tecnológicos de 
fabrico na fase de projecto. Deste modo, a utilização 
de modelos analíticos, implantáveis em programas 
gerais de análise de circuitos, torna-se desejável. 

Todavia, os modelos residentes nos simuladores de 
circuitos correntemente disponíveis referem-se a tran- 
sistores de junção bipolares (TJB) de baixa potência, 
como os que são implantados em circuitos integrados 
monolíticos. Tais modelos não só descrevem disposi- 
tivos estruturalmente diferentes dos TBPs, como tam- 
bém não incluem os efeitos dos principais mecanismos 
internos de alto nível que condicionam e limitam o 
funcionamento em sinais fortes. 

Uma parte importante do trabalho científico do 
Prof. Pedro Teixeira incidiu na caracterização física 
aprofundada de tais mecanismos internos, com vista à 
sua formulação, e no desenvolvimento de uma hierar- 
quia de modelos de circuito para o TBP de fácil inser- 
ção em simuladores de circuitos existentes. Além do 
mérito das contribuições originais na formulação e 
interpretação física dos fenómenos que caracterizam 
o funcionamento dos TBPs, deve salientar-se o modo 


particularmente bem sucedido como o autor assegurou 
a consistência interna das hipóteses que servem de base 
aos modelos para os diferentes níveis hierárquicos. 

Neste artigo, descrevem-se na Secção 2 as principais 
contribuições científicas do Prof. Pedro Teixeira para 
a modelação dos fenómenos físicos que condicionam o 
comportamento de TBPs. A caracterização da hierar- 
quia de modelos desenvolvidos e a sua implementação 
no simulador de circuitos SPICE2 constitui o objecto 
da Secção 3. Na Secção 4 demonstra-se como a estraté- 
gia de modelação de TBPs se pode generalizar e aplicar 
a dispositivos MOS (metal-óxido-semicondutor) de 
potência. Finalmente, apresentam-se na Secção 5 alguns 
exemplos de simulação de circuitos com TBPs, que 
evidenciam a eficiência e a precisão dos modelos e a 
sua utilidade no projecto de circuitos. 


2. MODELAÇÃO DOS FENÓMENOS FÍSICOS CA- 
RACTERÍSTICOS DOS DISPOSITIVOS BIPO- 
LARES DE POTÊNCIA 


2.1. Eficiência de Injecção para Níveis de Injecção 
Elevados 

Diversos autores apresentam a formulação das equa- 
ções de transporte no semicondutor como função explí- 
cita do mecanismo de transporte (deriva ou difusão) e 
do seu peso relativo, consoante o nível de injecção dos 
portadores de carga. Contudo, o autor desenvolveu a 
formulação das equações de transporte em termos dos 
quase-potenciais de Fermi, o que lhe permite caracteri- 
zar a operação de uma junção pn em regime de alto 
nível sem fazer intervir explicitamente as equações de 
transporte. 

Utilizando a aproximação regional de Schokley [1] 
c admitindo o quase-equilíbrio na região de transição 
da junção, verificou que o produto (pn) não depende 
da tensão aplicada à junção, Vs, mas da diferença dos 
quase-níveis de Fermi referentes a portadores maiori- 
tários entre as fronteiras de região de transição, APpu. 
Uma vez que só se verifica V; = APpx em baixo nível, 
propõe a utilização de 4Ppy como variável no modelo 
da junção, e relaciona APpy com a tensão exterior apli- 
cada, Va, e com as quedas de tensão ohmicas ao longo 
das regiões quase neutras n e p, Vo através de (ver 


fig 1) 


Fig. 1 — Junção pn polarizada directamente 


Fig. 2 — Variação de hp com J, 


Na qe Wiz + Va = AD — Vln [250/05 | —— 


— Vyln | Po Pr | + Vo (1) 


Esta formulação permitiu demonstrar a incorrecção 
de assumir como válidas as condições fronteiras de 
Fletcher [2] em alto nível, o que tem sido admitido 
para justificar (infundadamente) o estabelecimento de 
modelos empíricos; permitiu adicionalmente generali- 
zar as condições de Fletcher para qualquer nível de 
injecção. 

A eficiência de injecção da junção do emissor de 
um TBP e a sua dependência da excitação exterior (cor- 
rente que a percorre e tensão aplicada) determinam, 
em estruturas típicas, a lei de variação do ganho de 
ccrrente do TBP, hre, com a corrente de colector. 
O Prof. Pedro Teixeira propôs uma nova relação para 
a eficiência de injecção de uma junção que, para uma 
junção p*n se escreve Y = Jp/], ] = Jn + Jp, como [5] 


Y + (Lnp/Wn)a 


as ————— 2 
AM NA (Ln,/ Wa)a ve) 


p' (Xn) | (3) 


No 
P (Xn) | No + Pno 


nº (xp) n(Xp) Np | 


b = D,/D, = t,/&, e os demais símbolos retêm o seu 
significado usual. O parâmetro a reflecte o nível de 
injecção na região n (única onde poderá ocorrer injec- 
ção de alto nível), variando entre 1 (baixo) e 2 (alto 
nível). Atendendo a que ] = Jnp/(1 — 7), e que Jp = 
= q Dip nº (xp)/Lnp, pode obter-se a característica está- 
tica I(V) da junção exprimindo a e y (ou n' (xp) € 
p' (xn)) em termos de Va. 
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Para um TBP, o ganho de corrente em regime está- 
tico, hre, pode escrever-se 


Jp Y 


A = |—Y 


(5) 


her = 


o que permite obter uma lei de variação hre () = Jr) 
como a ilustrada na fig. 2. A lei de variação proposta 
permite verificar que Y passa uor um máximo, em 
regime de alto nível, o que está em harmonia com o 
facto de Y crescer (a y = c'*) com a e decrescer (aa = 
= c'*) com y. Permite ainda pôr em evidência o trans- 
porte adicional por deriva de lacunas na região n, em 
alto nível. 


2.2. Recombinação Directa no Emissor 


No funcionamento de um TBP, o nível de injecção 
condiciona o peso relativo de diversos mecanismos de 
recombinação sobre o comportamento terminal do dis- 
pesitivo. Assim, para níveis de injecção crescentes, 
tem-se sucessivamente (1) a recombinação SRH [4] 
(Shockley-Read-Hall), com um único nível de arma- 
dilhas, na região de transição da junção do emissor, (2) 
a recombinação por choque (Auger) [5] no corpo do 
emissor e (3) a recombinação SRH directa no volume 
da base e no epicolector [6]. O efeito da recombinação 
Auger no emissor faz-se sentir mesmo para a operação 
em baixo nível. Deste modo, tem de ser incluído na 
modelação da eficiência do emissor, uma vez que con- 
diciona o valor da corrente de portadores minoritários 
injectada no emissor. A recombinação Auger conduz 
a tempos médios de vida (ou a comprimentos de difu- 
são, Lpz) muito reduzidos. O autor propôs incluir a 
mcdelação do efeito Auger através de um factor cor- 
rectivo, bre, a ser aplicado à corrente minoritária do 
emissor, dado por 


di sea Toz a calo) a(o) 6 
ido lo qÃe WeNp(O) Ly: o, Lor (6) 


em que Qro é o integral da concentração na região 
neutra do emissor. Deste modo, hr, calculado com 
base no quociente dos números de Gummel do emis- 
sor e da base, vem reduzido de bar relativamente ao 
valor calculado sem este efeito. 


2.3. Modelo de Comando pela Carga Integral Gene- 
ralizada 


Num TJB de baixa potência, o integral da carga 
de base associada com o correspondente perfil de im- 
purezas, condiciona de forma dominante o comporta- 
mento do dispositivo. A percepção deste facto sem 
dúvida contribuiu para o estabelecimento da relação 
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de Gummel [7] e do modelo de Gummel e Poon [8], 
que procura descrever o efeito de alto nível através 
da variação das correntes terminais com a carga da 
base. Essa modelação é feita para a componente de 
transporte dominante de corrente entre o emissor e o 
colector; Gummel e Poon introduzem em seguida, de 
forma empírica, a corrente de base como perturbação 
das componentes de transporte. 


Todavia, o modelo de Gummel e Poon, implan- 
tado nas bibliotecas de simuladores de circuito como o 
SPICE 2 [9], é insuficiente para modelar correcta- 
mente o funcionamento de TBPs. Isto dá-se especial- 
mente para os regimes de saturação e quase-saturação, 
e deve-se à influência notável da carga integral asso- 
ciada à região do colector. O autor propôs assim uma 
generalização da relação de Gummel, em que a carga 
integral da base se encontra substituída pela soma das 
cargas integrais referentes às regiões que formam a 
junção, qsr ou qec, normalizadas em relação às respec- 
tivas cargas maioritárias na situação de equilíbrio, Qpso, 
Ouro € Quco. ÀS expressões propostas para as densidades 
de corrente longitudinais de electrões e de lacunas, nas 
fronteiras das regiões de transição, permitem determi- 
nar a densidade de corrente de base (ver Tabela 1 e 
fig. 3). 

A peneralização proposta pelo Prof. Pedro Tei- 
xeira, para além de conter, como caso particular, a 
formulação de Gummel, permite evidenciar (1) a não 
simetria do processo de transporte na base, sugerida 
por Gummel e que só é válida quando as cargas asso- 
ciadas a E e € são desprezáveis, (2) a expressão analí- 
tica das componentes da corrente de base, em alto nível, 
(3) a inclusão dos termos de recombinação, numa base 
não fenomenológica, (4) a extensão do efeito de Early 
em alto nível, (5) o mecanismo de retroacção através 
da modulação do volume de recombinação na base, e 
(6) a extensão da formulação do modelo à região do 
epicolector, o que torna possível descrever correcta- 
mente o desempenho do TBP em regimes de saturação. 


2.4. Desfocalização do Emissor 


A queda de tensão transversal na base, geralmente 
elevada na estrutura planar, típica dos transistores de 
potência para RF, provoca uma polarização decres- 
cente da junção BE, a partir do contacto de ligação da 
base. Deste modo, para um dada tensão exterior Ve, 
as correntes Is e Ic apresentam valores inferiores aos 
correspondentes à polarização uniforme da junção. 


EQUAÇÕES DO MODELO DE COMANDO PELA CARGA 
INTEGRAL GENERALIZADA 


Para a caracterização deste fenómeno, ao contrário 
de Hauser [10], Webster [11], ou Rey [12], o autor 


= MRE jd 
Ex qm b PE q 
J 
BS b+] “pn l b+1 
- | + ap E a  —— (] + a 
UBE b n; Ng pe 
' Q! 
. b + (b+1) q " [C1-Ap) s + AE — | 
E 


TABELA 1 - EQUAÇÕES DO MODELO DE COMANDO PELA CARGA INTEGRAL GENERALIZADA 
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Fig. 3 — Representação esquemática dum transistor n*pn—n 
unidimensional 


não estipula a geometria das linhas de corrente de base, 
pelo que os resultados obtidos são válidos para uma 
base não homogénea e incluem o efeito de modulação 
da condutibilidade da base, que ocorre em alto nível. 
No sentido de melhorar a simulação do regime dinã- 
mico, e de conferir um carácter bidimensional ao fun- 
cienamento na quase-saturação, apresentou um modelo 
resultante da divisão tranversal da base em duas 
secções. 


Modelo distribuído 


A partir do modelo de carga integral demonstrou-se 
que a iluminação do colector pela corrente de trans- 
porte em baixo nível é reduzida duma quantidade 
Fp (0) (factor de desfocalização), dado por 


Ep (0) bro senO cosO (7) 
O = e. eee 
" Iyp* 9 


em que o parâmetro O é função da corrente directa 
de base Iso, pesada pela corrente de referência Iro= 
=2Vr/Rspo, de acordo com: 


Igp 


OtgO = = k (8) 


Estando a corrente de transporte directa Irp expres- 
sa em termos da correspondente corrente do modelo 
de carga integral sem desfocalização (Rso=0), a mode- 
lação deste efeito é realizada através de um termo adi- 


cional àquele modelo. Em conformidade ,a sua implan- 
tação em programas de análise de circuitos que dispo- 
nham de modelo de carga integral é trivial [13]. 

Para correntes elevadas, a concentração de electrões 
iniectados pelo emissor na base torna-se comparável à 
concentração de lacunas em equilíbrio. Por conse- 
guinte, a resistência tranversal da base descresce para 
valores crescentes de Isp (modulação da condutibili- 
dade da base). Trata-se de um efeito que contraria e 
atenua a desfocalização de emissor, em alto nível. 

O Prof. Pedro Teixeira considerou a carga da base 
normalizada, qs, constante segundo a direcção trans- 
versal da base e definida pelo modelo de carga integral 
sem desfocalização; deste modo, a introdução da modu- 
lação da condutibilidade da base num transistor com 
desfocalização reduz-se à substituição em (8) da cor- 
rente de referência em baixo nível Iro por Ia=2Vr/Rs 


onde Rg é a resistência da base em alto nível, definida 
como: 
Reo 
Ná = (9) 
qB 


Refira-se que a técnica de modelação utilizada per- 
mite uma generalização fácil para outras regiões de 
funcionamento, nomeadamente a zona inversa (desfo- 
calização de colector) e a de saturação. 

Para verificar a hipótese que conduziu a (9), o 
autor recorreu à simulação dos efeitos bidimensionais 
poi um modelo em cadeia de elementos concentrados. 
Para tanto, considerou-se o dispositivo dividido num 
número elevado de secções, ao longo da dimensão trans- 
versal (fig. 4). Na figura 5 apresenta-se a comparação 
entre as características normalizadas (Va=V/Vr, In= 
=| Rso/Vr, Rn=R/Rso), Ic(Vbe) obtidas para o tran- 
sistor 2N2222 a partir do modelo distribuído (7-9) e 
do modelo de elementos concentrados com 100 secções. 


Modelo de elementos concentrados — 2 secções 


Na simulação do regime dinâmico, a capacidade 
do programa de análise de circuitos e o tempo de com- 
putação vão limitar o número de secções a utilizar pela 
representação de elementos concentrados. Em confor- 
midade, o Prof. Pedro Teixeira estabeleceu um modelo 


Fig. 4 — Cadeia de N secções — com modulação de condutibilidade da base 
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— N = 100 


o mod. distr. 
(7 ='5) 


VBEn 
25 30 35. 40 


Fig. 5 — Características normalizadas 1 (Vgp) e Ia(Vpg) com 
modulação da condutibilidade da base 


Fig. 6a — Partição em duas secções transversais 


Fig. 6b — Esquema de princípio 


de duas secções que verificou conduzir a boas aproxi- 
mações da solução exacta (fig. 6). Por comparação com 
a resposta em regime transitório, dum modelo distri- 
buído aproximado por um número elevado de secções, 
definiu-se um critério para a definição do parâmetro 
de partição a (a = 0.15 — 0.30). 


2.5. Efeitos do Epicolector 


Os efeitos associados à região do colector são geral- 
mente dominantes na operação em alto nível dos TBPs. 
No seu funcionamento como elemento do circuito, veri- 
fica-se a existência de uma zona intermédia de quase- 
-saturação (ou saturação parcial) nas características 
estáticas de saída, não prevista pelo modelo de carga 
integral. 

Há dois efeitos interactuantes em alto nível que 
influenciam o comportamento do TJB na região de 
quase-saturação: a penetração da região efectiva da 
base para o interior do epicolector n e a desfocaliza- 
ção do emissor já referida atrás. As características vin- 
cadamente bidimensionais do fenómeno de desfocaliza- 
ção de emissor arrastam consigo a necessidade de con- 
duzir também o estudo do fenómeno de alargamento 
de base numa perspectiva bidimensional, em especial 
para a correcta modelação de regimes transitórios. 
O modelo de duas secções tranversais proposto pelo 
autor revelou-se, todavia, suficientemente preciso para 
a caracterização em vista. No entanto, a modelação por 
uma única secção conduz também a resultados com 
precisão suficiente na simulação e é normalmente adop- 
tada por razões de eficiência de computação. 

Deste modo, o alargamento da região de base efec- 
tiva, provocada pela injecção de portadores minoritá- 
rios no epicolector n, é o fenómeno dominante da 
quase-saturação. Ele produz essencialmente três efeitos 
ra operação do dispositivo: a modulação da resistên- 
cia de corpo do epicolector, a redução do ganho de cor- 
rente hre e, para o regime dinâmico, os efeitos relacio- 
nados com o armazenamento de cargas no epicolector. 


Resistência do epicolector 


A partir das equações de transporte no epicolector, 
suposto dividido em duas regiões, uma de injecção 
(o—x;) e outra de condução (x;— Wc), onde x; é defi- 
nido a partir da condição de equilíbrio p(x;) =n//Npe, 
o autor obteve a importante relação: 


VcB = [Rco—Kc E Vaci/Vr) — | = Vec; (10) 


onde Re=2Vr (n;/Npc”) e Rco é a resistência ohmica 
do corpo do epicolector. 

Em termos de implantação num programa de aná- 
lise de circuitos, a expressão (10) corresponde a uma 
malha com um TJB elementar (Vecj), uma resistência 
ohmica Rco e a uma resistência não linear de colector 
simulada por um gerador comandado por Vec;. 


Ganho de corrente hrr 


Dois mecanismos são responsáveis pelo decréscimo 
do ganho de corrente na quase-saturação: o aumento 
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da carga de portadores maioritários na base, devido ao 
alargamento, e a recombinação no colector. 

Sendo o primeiro mecanismo responsável por um 
decréscimo Alc da corrente de colector e o segundo pela 
corrente inversa de base Isr, o ganho de corrente hre 
é dado por: 


| —(Alc/Lro) | 
(11) 


Bi) = 
da “Ep | 1 + (Tor/Tap) 


Estes efeitos foram introduzidos no modelo global 
de uma forma conspícua. Para tanto, consideraram-se 
correntes aditivas ou subtractivas que incidem na cor- 
rente total de colector, relativamente ao valor previsto 
pelo modelo de transporte em baixo nível (lo, 
lap) [14]. 


Regime dinâmico 


Para a descrição do regime dinâmico na quase- 
saturação, incluiu-se um termo adicional na corrente 
de base, que exprime a variação no tempo da carga 
armazenada no epicolector Qpc(t). Há portanto que 
introduzir no modelo um gerador adicional entre base 
e colector, que fornece a corrente: 


dOpclt) 


dt (12) 


Para a caracterização de Qpc(t) foi adoptada a apro- 
ximação de comando pela carga. Foi também introdu- 
zida uma aproximação linear para sinais fracos a partir 
da relação: 


GQ — — Csveb + Trele (13) 


onde os parâmetros Tre e Cs são dependentes da pola- 
rização (le, Vce). 


3. MODELO PARA O TBP 


5.1 Hierarquia de modelos 


A partir da caracterização com base física dos efei- 
tos de alto nível referida na secção anterior, o Prof. 
Pedro Teixeira desenvolveu um modelo analítico para 
o TBP para implantação no programa para simulação 
de circuitos SPICE 2 ou similares. O modelo foi inse- 
rido sob a forma de um circuito equivalente constituído 
por diversos subcircuitos separáveis, estando cada um 
er: princípio ligado à caracterização de um efeito de 
alto nível. 

Tornou-se assim possível conceber uma hierarquia 
de modelos, de complexidade crescente, com base na 
sebreposição (da incidência parcial) dos mecanismos 
físicos mais relevantes no funcionamento em alto nível. 
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EFEITOS 


Acção elementar do transistor 
Armazenamento de cargas na base 


| Efeitos parasitas laterais 


| Nível I 


Injecção em alto nível (ambipolar) na base 
| Armazenamento de cargas no emissor e colector 
Componentes de recombinação explicitadas 


Modulação do temo de trânsito 


Nivel II 

Efeitos estáticos e dinâmicos do colector na 
saturação parcial e ccmmleta 

Quedas de tensão internas 

Desfocalização de emissor em alto nível 
Disrunção 


Simulação bidimensional por secções transversais 


TABELA 2 


Estes mecanismos encontram-se ligados entre si, e a 
sua ocorrência depende das condições exteriores de 
polarização. Os níveis hierárquicos adoptados (I, II 
e 111) correspondem a um certo compromisso entre a 
verosimilhança do modelo, a complexidade da implan- 
tação e o tempo de cálculo associado à simulação em 
computador. A repartição dos diferentes efeitos por 
nível de modelação está descrita na Tabela 2. 

As expressões mais relevantes para a caracterização 
dcs referidos efeitos estão, por sua vez, sumarizadas 
na Tabela 3. Esta tabela inclui também as relações refe- 
rentes ao transistor elementar (nível 1), o qual funciona 
como elemento de ligação na modelação parcial daque- 
les efeitos. (Considera-se a operação na zona directa). 


Uma estatística da hierarquia proposta, referente 
a um circuito padrão (inversor com carga RL), é apre- 
sentada na Tabela 4. No sentido de melhorar a eficiên- 
cia da simulação, o Prof. Pedro Teixeira propõe ainda 
a adopção de um macromodelo funcional [13] que 
contempla, em particular, a operação na saturação. 


3.2. Parâmetros de limiar 


No sentido de obter resultados consistentes, pode 
electuar-se uma selecção prévia do nível do modelo a 
adoptar através da aplicação de regras expeditas de 
dimensionamento e projecto. Este processo de escolha 
pode porém sistematizar-se a partir da introdução de 


TABELA 3 
E Nº nos 
Nivel 
internos 


AOC 
TABELA 4 
parâmetros de limiar exclusivamente dependentes da 
estrutura do dispositivo. Os valores destes parâmetros, 
determinados com base na caracterização dos efeitos a 
que dizem respeito, fixam um limite inferior (no domi- 
nio das variáveis eléctricas aos terminais) para uma 
incidência notável do respectivo efeito quando conside- 
rado em separado. 

A determinação da hierarquia requerida para o 
modelo completo do TBP é imediata no caso da separa- 
ção pronunciada dos valores dos parâmetros de limiar. 
Por exemplo, se Ians<Ic<los<Ice, as características 


Nº derelaçoes 


nao lineares 


DISRUPÇÃO DIRECTA (X) 
|KxVa<Vcgs Vi | 


Tempo médio 
deCrU(s) 


É 


aos terminais nestas condições de funcionamento vão 
reflectir, de forma predominante, os efeitos de injecção 
de alto nível na base, nomeadamente quanto à depen- 
dência de Irn(Vse) e hre(Vpr). Pode, por conseguinte, 
empregar-se o modelo de nível II, na simulação do seu 
comportamento num circuito de aplicação em que se 
verifiquem a priori aquelas desigualdades. 

No entanto, no funcionamento dum dispositivo real, 
vai ocorrer a interacção dos diversos efeitos, a qual 
assume fundamentalmente dois aspectos distintos: fora 
da saturação ou da quase-saturação, há sobreposição 
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de CE e ANB, que se pode caracterizar de forma unidi- 
mensional por uma relação do tipo Ice=2plro para a 
componente directa de transporte. A partir do limiar 
da entrada na quase-saturação (na periferia da base) 
ocorre guase-saturação parcial, que requer uma carac- 
terização bidimensional. A condição de ocorrência da 
quase-saturação parcial foi verificada como sendo 
[c>Ios=Fp los, em que Fp é o factor de desfocaliza- 
ção e los a corrente limiar de quase-saturação. A região 
transversal quase-saturada passa por um valor máximo, 
tendendo assintoticamente para zero. O estudo numé- 
rico do modelo bidimensional permite concluir que 
aquele máximo pode ser igual ou maior que a dimensão 
transversal do emissor. No caso de yo 1 (quase-satu- 
ração completa), é válida a descrição unidimensional. 
Para valores crescentes de Ic tem-se portanto a seguinte 
sucessão: 


— [e<Tos — zona activa directa; 

— [cos zona de quase-saturação parcial 
(ro<]); 

— eventualmente (caso yo 1), uma zona intermé- 
dia de saturação completa; 


— uma região assintótica de quase-saturação par 
cial (yo ”* 0). 


Concluído assim que a quase-saturação do disposi- 
tivo é parcial ao longo de uma faixa extensa de valores 
(1,V), foi possível estabelecer, para o nível III, a hie- 
rerquia de modelos da Tabela 5 (não considerando a 
disrupção). 


O modelo do nível 111.3 pode incluir apenas OS, 
uma vez que se trata do efeito dominante para aquela 
faixa de correntes. Deve no entanto referir-se que, em 
muitos casos, o modelo unidimensional EC+ QS+ ANB 
assegura uma caracterização suficientemente precisa, 
em especial quando a estrutura do dispositivo conduz 
a valores pouco elevados da resistência transversal da 
base em equilíbrio. 


3.3. Caracterização dos modelos. Determinação dos 
parâmetros 


Os parâmetros de implantação dos modelos analíti- 
cus apresentados atrás encontram-se ligados à descrição 


nivel TJ.) CE + ANB 


nivel III.2 | CE +ANB +QS 


nivel III.3 CE + ANB + QS 


física de um certo número de efeitos. À sua dis 
nação pode ser efectuada: 


— com base na caracterização fina, física e dimen- 
sional da estrutura; 

— a partir de um programa consistente de medidas 
(estáticas e dinâmicas) aos terminais do dispo- 
sitivo real encapsulado ou directamente sobre a 
pastilha. 


A incidência dum ou doutro modo vai depender do 
critério da utilização dos modelos: em circuitos de 
aplicação, no projecto de dispositivos ou no controlo 
de processos. Nos trabalhos do Prof. Pedro Teixeira, a 
caracterização foi orientada no sentido da identifica- 
ção completa dos parâmetros através dum programa 
de medidas aos terminais. O autor apresentou, no 
entanto, exemplos que verificam a concordância entre 
os valores dos parâmetros resultantes das medidas com 
os obtidos a partir duma descrição dimensional e física 
aproximada. Apresentam-se, na Tabela 6, com indica- 
ção dos níveis de modelação, os valores numéricos para 
os 39 parâmetros do modelo para três transistores de 
aplicação corrente: 2N918 (sinais fracos, alta fre- 
quência), 2N3632 (média potência, alta frequência) e 
BUW34 (Darlington de alta potência para comutação), 


4. EXTENSÃO À MODELAÇÃO DE TRANSISTO- 
RES MOS DE POTÊNCIA 


O processo de comando da corrente nos transistores 
MOS é diverso dos bipolares; no entanto, ambos apre- 
sentam uma região de baixa condutividade (epidre- 
no/epicolector) em série com a região activa, por forma 
a suportar altas tensões (fig. 7). Deste modo, é possível 
caracterizar o funcionamento dos TECMOS de potên- 
cia em sinais fortes por uma técnica idêntica à utilizada 
para os bipolares nos parágrafos anteriores [14]: mode- 
lação, em blocos independentes, dos principais efeitos 
que dominam o comportamento aos terminais, nomea- 
damente o comando da corrente pelo terminal da porta 
e a resistência não linear do epidreno [15]. 


O paralelismo existente entre a técnica de modela- 
ção utilizada pelo Prof. Pedro Teixeira para os TBP 
e a seguida pelo seu colaborador J. Costa Freire para 


unidimensional 


duas (ou mais) secções 


unidimensional 


TABELA 5 
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TABELA 6 


Designação do grupo | nível do modelo 


tema 
Ea! 
ld 


modelo básico de transporte 


tempos de trânsito intrínsecos 


tempos de transito laterais 


a caracterização da operação dos VVMOS ()) em sinais onde é utilizada a mesma topologia para as malhas de 
fortes [16] é realçada pelo esquema básico de ambos entrada e saída, embora com diferentes significados 
os transistores de potência representado na fig. 8, para as respectivas variáveis. 


(1) VVMOS: Transistor MOS de canal Vertical e de estrutura com cunha em Ng. 7. 
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Fig. 7 — Estruturas típicas (a) de um transistor VVMOS e 
(b) de um TBP 


Refira-se que, embora a implementação de ambos 
os modelos seja similar em termos de circuito básico 
equivalente, o modelo do TBP é consideravelmente 
mais elaborado, já que a condução é essencialmente 
ambipolar, envolvendo fenómenos de alto nível e de 
armazenamento de portadores minoritários, como se 
referiu anteriormente. Este facto é salientado na fig. 9, 


|C'| 

E RS | 
ad E | 

l bp 

É Mis 

L epl 
“pl o 

| | 

| | 

Ê l 


, = mm o OS SO 


onde se representa de uma forma simplificada os 
circuitos de implantação no programa SPICE 2 dos 
modelos desenvolvidos para estes dois tipos de transis- 
tores de potência. 

Para o TBP introduzem-se adicionalmente ao mode- 
lo de carga integral residente duas resistências não 
lineares associadas aos fenómenos de desfocalização 
de emissor (DFE) e de modelação da resistência de 
colector (RC), um gerador de corrente (Mic) que intro- 
duz o efeito de variação do ganho hr: na quase-satura- 
ção e uma capacidade não linear (Qpc) referente ao 
armazenamento de portadores minoritários no epicolec- 
tor. Para o VVMOS é apenas adicionada ao modelo de 
TECMOS de canal curto, uma resistência não linear 
(RD) que tem em conta a dispersão da corrente no 
epidreno e a sua dependência com a configuração do 
canal junto ao epidreno, e a dependência da mobili- 
dade, nesta última região, com o campo eléctrico. Para 
a caracterização do epidreno foi utilizada uma técnica 
de partição regional idêntica à seguida pelo Prof. Pedro 
Teixeira para o TBP (secção 2.5). 


5. EXEMPLOS 


Referem-se seguidamente dois exemplos de simula- 
ção apresentados pelo Prof. Pedro Teixeira para provar 
a validade e precisão do modelo por si proposto para 
o TBP, 


5.1. Amplificador de RF de classe €C para HF 


Este primeiro exemplo corresponde a um caso de 
diagnóstico: o amplificador de RF de classe C para 
[=7MbHz, cujo esquema preliminar se apresenta na 
(fig. 10a) foi correctamente projectado com base nos 
dados de catálogo. 


h 


reltV ro Vos ig) 


Relio: Vepi! 


Fig. 8 — Modelo de circuito elementar para TBP e VVMOS e definição de parâmetros 
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Fig. 9 — Modelos para o TBP e para o VVMOS implantados 


Na construção de diversos protótipos foi possível 
observar um rendimento de conversão baixo (1 =30— 
— 40%) e a destruição do TBP por disrupção primária 
da junção base-colector (vce, > Vc, ). Apresentam-se 
na fig. 11 os oscilogramas de ve(t) e vo(t) para aquele 
regime de funcionamento precário. 


Vett) 


Volt) 


0+13.6V 


| Ls 
2.2 nF = 100nF El 
SuH E L 
| St 
O.luF 0.5uH | 
co | 


ve ft) a 
280 pF 


| 


o— o—— o— ——— O um ] 


no simulador SPICE 2 


A simulação do circuito da fig. 10, recorrendo ao 
modelo de nível III, forneceu para vc(t) e vo(t) a evolu- 
ção que se reproduz na fig. 12 (regime permanente). 
A variação no tempo de vc(t) é caracterizada por uma 
componente de alta frequência ([=66MHz), que ocorre 
por salvas, com uma frequência de repetição igual à do 
sinal de entrada. O valor de vc(t) pode, nestas circuns- 
tâncias, exceder Vce,.. indicado pelo fabricante. 


1 


Os resultados anteriores permitiram concluir não 


o a a dr 


Fig. 11 — Oscilogramas de volt) e volt) — amplificador 7 MHz classe C 
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Fig. 12 — Resultados da simulação do amplificador de 7MHz, 


se estar em presença de uma instabilidade propriamente 
dita, mas sim de um caso particular de excitação da 
malha de adaptação de saída ocorrendo no período de 
corte do TBP. Durante este período, a corrente de ali- 
mentação Ic (através de L;) vai excitar o circuito de 
saída, cuja frequência de ressonância é, nesta configu- 
ração, determinada por L; e pela capacidade total de 
saída do transistor Cr, a qual inclui as capacidades 
parasitas de montagem, de acordo com o esquema da 
fis. 10b). A forma de onda de vc(t) é, por conseguinte, 
dependente do tipo de malha de adaptação de saída 
empregue no amplificador. Os seus efeitos podem, do 
ponto de vista do valor máximo de vce(t), ser minimi- 
zados aumentado Cr através da inclusão de uma capa- 
cidade adicional entre o colector e a massa de forma a 
que o novo valor de Cr sintonize L; para uma frequên- 
cia próxima da de trabalho. Dada a configuração da 
malha, as alterações a efectuar em L; e C; para reajus- 
tamento da adaptação à carga são mínimas. 

Foi assinalado ainda pelo Prof. Pedro Teixeira o 
facto de os resultados da simulação revelarem não 
serem iguais as amplitudes máximas de vc(t) em duas 
salvas consecutivas, o que sugere uma modulação de 
amplitude adicional. Este efeito não ocorre na simula- 
ção com os modelos de nível | ou II, pelo que pode 
estar relacionado com qualquer efeito (por exemplo, 
paramétrico) de colector, considerando o autor ser este 
um assunto em aberto. 


5.2. Amplificador (conversor) de classe E 


Apresenta-se na fig. 13 o esquema de ligações dum 
amplificador comutado de classe E. Neste tipo de 
amplificador, a malha de saída é dimensionada de for- 
ma que uma das variáveis vclt) ou ic(t) seja nula 
quando a outra o não fôr. Consegue-se, deste modo, 


a4 


4 


800 900 


1000 


classe € (modelo de nível 111) 


atingir rendimentos de conversão muito elevados, tipi- 
camente da ordem dos 90%. 


Num caso real, o dimensionamento da malha e o 
desempenho da montagem dependem fortemente das 
características estáticas e dinâmicas do TBP, em espe- 
cial das associadas à saturação completa ou parcial, 
uma vez que o dispositivo funciona como interruptor. 
E por conseguinte vantajoso o emprego do modelo de 
nível III na simulação. Os valores dos elementos da 
malha de saída foram calculados de acordo com o 
método proposto por Raab [17] que considera, por 
simplicidade, o TBP no circuito como um interruptor 
ideal. Procedeu-se seguidamente a uma optimização 
daqueles valores, através de simulações repetidas, no 
sentido de assegurar um rendimento de conversão ele- 
vedo. Apresenta-se, na Tabela 7, a comparação entre 
os valores calculados, optimizados na simulação, e os 
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Fig. 13 — Amplificador comutado de classe E 


TABELA 7 


Variável 


1 
1 
1.85 


4 
0 
356 

5 


7 
1 


Calculado 


Simulação Experimental 


0 


valores experimentais de algumas das grandezas e de 
componentes associados ao circuito da fig. 13. A 
figura 14a) mostra os oscilogramas de vclt) e ic(t); 
a fig. 14b) mostra os resultados da simulação do cir- 
cuito (versão final) para aquelas grandezas. 


a) Oscilograma de ic(t) e de veglt) 


b) Simulação com os modelos de níveis [ e HI 


6. CONCLUSÕES 


Apresentou-se neste artigo uma colecção das prin- 
cipais contribuições originais do Prof. Pedro Teixeira 
para a modelação e simulação de dispositivos activos 
semicondutores, tendo como destino o projecto de cir- 
cuitos electrónicos. Para os TBPs, evidenciou-se a mode- 
lação da eficiência de injecção da junção do emissor, 
da recombinação directa no volume da região quase- 
neutra do emissor, dos efeitos de desfocalização no 
emissor e de injecção em alto nível na região do epico- 
lector e do regime de operação em saturação parcial. 
Salientou-se o modelo de comando pela carga integral 
generalizada proposto pelo autor, e que contém como 
caso particular o modelo de Gummel e Poon. 

Do ponto de vista de simulação, apresentou-se uma 
hierarquia de modelos proposta pelo Prof. Pedro Tei- 
xeira para apoiar o projecto de circuitos de potência, 
definindo-se o domínio de validade de cada modelo em 


função das variáveis eléctricas terminais e a forma de 
determinação dos seus parâmetros. 
Os exemplos de simulação descritos procuram 
salientar a utilidade dos modelos desenvolvidos para o 
0.4 = eus 


Fig. 14 — Resultados obtidos para o amplificador comutado 
de classe E 
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projecto de circuitos, e a indispensabilidade de tomar 
em consideração efeitos tidos como de segunda ordem 
em modelos anteriores, não só para a optimização do 
desempenho dos TBPs como também para garantir o 
seu funcionamento seguro. 

Finalmente, procurou-se dar relevo ao facto de que 
a estratégia de modelação que o autor empreendeu 
pode aplicar-se a outros dispositivos activos, nomeada- 


mente 


a transistores VVMOS de potência, 
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RESUMO 


As evoluções no conhecimento e no valor real da energia 
estão a convergir no sentido de tornar o uso das energias 
alternativas concorrencial com as soluções clássicas para 
fornecer energia eléctrica a dreas rurais isoladas. Para si- 
tuar as soluções alternativas perante as soluções clássicas 
estabelecem-se alguns valores que actualmente caracteri- 
zam estas últimas. Depois de estimar as utilizações anuais 
que se podem esperar das fontes eólicas e das fontes foto- 
voltaicas, confronta-se a qualidade da energia “bruta” ob- 
tida nas novas fontes com as exigências dos utilizadores. 
Disto resultam os requisitos genéricos para as interfaces, A 
comunicação termina com uma sugestão sobre critérios de 
avaliação do interesse dos investimentos. 


1. DEFINIÇÕES 
1.1 — Fontes de energia alternativas 


Por “fontes de energia alternativas” entendem-se apro- 
priações no fluir das energias naturais usando processos 
tecnológicos não clássicos. Em vez de “fontes de energia 
alternativas” seria mais apropriado escrever “novas formas 
de inserção no fluir natural da energia”, Contudo a expres- 
são “energias alternativas” está a generalizar-se universal- 
mente; deve esclarecer-se que a expressão exclui a energia 
nuclear. 

Entre as muitas fontes de energia alternativas que, hoje, 
abundam na literatura técnica — embora abundem menos 
as realizações concretas — será dada atenção particular, 
nesta comunicação, aos aerogeradores e às células fotovol- 
taicas. Afigura-se que são estas as fontes alternativas que, 
actualmente, merecem mais atenção como futuras compe- 
tidoras das soluções clássicas, quando se trata de abastecer 
com energia eléctrica áreas rurais isoladas. As próprias cen- 
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trais mini e microhidroeléctricas parecem só ter interesse, 
para este fim, em situações muito particulares. Dai o pouco 
destaque que, nesta comunicação se dá à energia de origem 
hidrica embora muito do que adiante se escreve lhe seja 
aplicável. 

Por “redes eléctricas rurais isoladas” entendem-se as que 
servem centros de consumo em áreas rurais distantes das 
redes eléctricas de distribuição regionais. As potências de 
ponta dos diagramas de carga nestas redes isoladas situam- 
-se entre um número digito e poucas centenas de qui- 
lowatts. 

Excluem-se pois desta designação as redes que servem 
aglomerados populacionais com muitas centenas ou mais 
de um milhar de habitantes; excluem-se, também, as cargas 
isoladas de pequenissima potência tais como algumas esta- 
ções repetidoras de sinais de telecomunicação, postos tele- 
fónicos ou postos meteorológicos isolados, sinais para a 
navegação e análogas. 

Por “interfaces” entendem-se os conjuntos de equipa- 
mentos eléctricos e electrónicos montados entre a fonte 
primária (veio da turbina eólica ou terminais da célula fo- 
tovoltaica, por exemplo) e a rede eléctrica. 

Nos centros de consumo em áreas rurais isoladas diver- 
sas circunstâncias concorrem para que a energia eléctrica 
seja usada apenas para fins em que, sem incorrer em gran- 
des incomodidades, não é possível substituí-la por outras 
formas de energia. Entre essas circunstâncias avultam duas: 
a) qualquer que seja a fonte primária usada, a unidade de 
energia eléctrica é muito cara; b) os orçamentos dos con- 
sumidores são modestos e não se deverá esperar que reser- 
vem parcela importante desses parcos orçamentos para gas- 
tos energéticos. Toda uma evolução se processou ao longo 
de séculos que explica que a área seja “rural” e “isolada”, é 
algo ingênuo pensar que o fornecimento de energia eléc- 
trica poderá, só por si, inflectir essa evolução. 

Acresce que as produções locais de energia eléctrica, com 
instalações de pequena potência e baixa utilização, não po- 
dem garantir qualidade da energia tão elevada como o pode 
fazer uma grande rede. É certo que não seria razoável espe- 
rar que, quanto à energia, se conseguissem, nessas áreas 
isoladas, altos padrões de qualidade que não se conseguem 
nem nos serviços de saúde, nem na educação, nem na cul- 
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tura, nem nos transportes. Todavia o certo é que a irregula- 
ridade no fornecimento da energia eléctrica, as oscilações 
na tensão € na frequência e outros defeitos, inevitáveis nas 
produções locais, não atraem os consumos, 

Contudo deve ponderar-se que, se a regulação de tensão 
e a taxa de harmônicas admissíveis para muitos dos recep- 
tores instalados na rede rural isolada devem ser análogas às 
que se obtêm nas grandes redes, já não se afigura necessária 
nem a apertada regulação de frequência característica da 
rede europeia interligada (+ 0,1] Hz) nem a elevada garan- 
tia de serviço que nela se verifica. 

A fixação dos padrões de qualidade minima a exigir à 
energia eléctrica fornecida pelas redes rurais isoladas 
aguarda ponderado e minucioso trabalho de investigação, 
sem o qual o desenvolvimento das várias soluções tecnoló- 
gicas anda à deriva, movido apenas pelos interesses de cada 
proponente ou pelas divulgações que deslumbram e enga- 
nam os financiadores. Neste particular o visar qualidade 
desnecessária pode ser tão prejudicial como aceitar quali- 
dade insuficiente. 

O conhecimento avança, o que tem como consequência 
baixar o custo real das energias alternativas; o valor real da 
energia é cada vez maior. Ás evoluções do conhecimento e 
do valor da energia convergem no sentido de tornar o uso 
das energias alternativas concorrencial com as soluções 
clássicas, cada vez em maior número de situações, 

O duplo ataque do problema, simultaneamente pela face 
tecnológica e pela face económica é por isso, o objecto 
desta comunicação. 


1.2 — Alguns números 


Apenas para situar a ordem de grandeza dos valores da 
potência, do consumo anual e da utilização que se podem, 
genéricamente, esperar nas redes rurais isoladas, considere- 
-se o exemplo seguinte. 

Numa área rural isolada concentram-se vinte potenciais 
consumidores domésticos. Um inquérito prévio leva a ad- 
mitir que a potência de ponta total será condicionada pelas 
seguintes potências parcelares, por consumidor: 


5 lâmpadas incandescentes de 60 W 300 W 
| televisor a cores (importante 

para as acções educativas) 140 W 

| frigorífico 260 W 

Total 700 W 


A potência de ponta, nos vinte consumidores, será, assim 
Pm = 14 kW. 

Pode admitir-se que o consumo anual, por consumidor, 
totalizará cerca de 2000 kWh; para este consumo concor- 
rem mais alguns electrodomésticos que se supõe, não esta- 
rem ligados na hora de ponta, tais como ferro de engomar, 
aspirador, torradeira, rádio-receptor e outros. A eventual 
bombagem de água para usos domésticos poderá também 
concorrer para o consumo de energia eléctrica. 

O consumo, nos vinte consumidores, será, então, 
Wa= 40000 kWh / ano. A utilização da potência de 
ponta será, então, h, = W, / Pm = 2860 horas por ano. 

Pode concluir-se que as redes rurais isoladas se caracteri- 
zam por utilizações mediocres, crescimento de consumo 
lento e imprevisível e exigências modestas quanto à quali- 
dade da energia eléctrica. 
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2. SOLUÇÕES CLÁSSICAS 
2.1 — Ligação à rede eléctrica regional 


As soluções clássicas para o fornecimento de energia a 
redes rurais isoladas são a construção de um ramal até à 
rede regional a média tensão ou a instalação de grupos 
diesel-eléctricos. Outras soluções têm sido usadas mas cor- 
respondem a decisões singulares explicadas (nem sempre...) 
por condicionamentos pouco frequentes. 

A ligação à rede de distribuição regional estabelecendo 
um ramal a média tensão é, na generalidade das situações, a 
melhor solução tecnológica. Contudo o custo da unidade 
de energia pode anular o interesse da solução. 

O custo C, de cada quilowatt-hora fornecido na extre- 
midade receptora do ramal pode ser avaliado pela expres- 
são 


1 4 W, 
É = +E.(1 + ) [escudos / kWh] (1) 
a u 
em que 

I, — Investimento total no ramal e terminais 
[escudos] 

| encargos de capital 
[escudo . escudo !.ano !] 

] — preço de compra da energia na emissão 
[escudo . kWh |] 

W, -— consumo anual [kWh . ano !] 

A Wa — energia perdida no ramal e terminais 


[kWh . ano !] 


O preço de compra da energia na emissão será regido, mais 
ou menos explicitamente, pela lógica das tarifas binómias. 
Todavia, para as redes rurais isoladas, será o primeiro 
termo do segundo membro da expressão | que avantajará 
no custo da energia. 

De facto as médias tensões permitem transmitir a distân- 
cias de alguns, ou de algumas dezenas de quilómetros, po- 
tências de centenas ou milhares de quilowatts (Sobre estas 
possibilidades é esclarecedora a sistematização feita no dia- 
grama que se reproduz em anexo e que é extraido da refe- 
rência [1]). Pode assim concluir-se que o investimento 1, 
feito para se ligar a rede geral à distribuição regional é 
sub-aproveitado. verificando-se 


ss 
a o Ss ff 4 
W 


a a 


O ramal trabalha praticamente em vazio sendo minimas 
as perdas na transmissão À W,. Tudo se passa como se o 
custo C, dependesse fundamentalmente do número total W, 
de unidades de energia veiculadas pelo ramal, sendo secun- 
dária a influência do valor da ponta P,, que afecta f. 

Se, para alimentar a rede rural exemplificada na alínea 
1.2, foi construido um ramal trifásico a 15 kV com 15 km de 
extensão, a expressão | conduz a custo C| = 70 escu- 
dos. kWh! , Para obter este custo admitem-se i = 0,25 
e f = Bescudos. kWh! Um preço f baixo da 
energia na rede de distribuição regional, não é condição 
suficiente para se poder oferecer energia eléctrica barata às 
comunidades rurais isoladas. 


É evidente que o aumento do consumo W, pode mudar 
esta situação e tal possibilidade deve ser examinada com 
cuidado. Todavia nas áreas rurais isoladas o crescimento 
do consumo é, em geral, imprevisível e provavelmente 
muito lento, como já se notou. 


2.2 — Central diesel-elécirica 


Entre as máquinas térmicas que garantem razoável con- 
tinuidade no serviço, o motor diesel é a mais barata por 
quilowatt instalado, isto para as pequenas potências ade- 
quadas às redes rurais isoladas. Explica-se assim a sua lar- 
guissima difusão que é explicada, também, por se dispor, 
em qualquer local habitado, de rede de distribuição de 
combustivel e de assistência competente: a rede e a assistên- 
cia que servem os motores dos transportes rodoviários de 
mercadorias. 

No Quadro | indicam-se custos unitários actuais (1983) 
de centrais diesel-eléctricas de pequenas potência. Os custos 
incluem os equipamentos auxiliares e a construção civil. 


QUADRO 1 — Custos de centrais diesel-eléctricas (1983) 


sda Velocidade 
Potência grs 
angular Custo unitário 
instalada 

dos grupos 

P; n l, 
[kVA] [rot.min!] | [10º escudos /kW] 
E) [500 70 
| x 125 1500 45 
| x 625 [500 40 
| x 625 LOD0 TO 
| x 1250 [500 40 
| x 1250 LODO 70 


Admita-se que, para alimentar a área rural isolada que se 
descreve na alínea 1.2 o investimento na central era 
|| = 80.10º esc/kW, que a utilização era h, = 2860 hora 
ano !; admita-se ainda que por cada quilowatt-hora forne- 
cido ao consumidor se gastavam 0,3 kg, de gasóleo ao preço 
de 30 esc/kg, (f.o.b. Roterdam, | dólar USA < > 120 
escudos, 1983). Então o custo unitário da energia será 


C, = l, ' ha + 9 = [6 escudos/kWh (2) 
isto se for | = 25% ao ano a taxa de encargos com o inves- 
timento. 


O exame pormenorizado da economia das centrais diesel- 
eléctricas, em áreas rurais isoladas, não é o objecto desta 
comunicação. Todavia para melhor situar as fontes de 
energia clássicas perante as fontes de energia alternativas 
convém referir ainda o seguinte: 


a) Para garantir razoável continuidade no fornecimento 
deverá dividir-se a potência instalada por dois ou três gru- 
pos. Isto elevará o investimento para cerca de 100.000 es- 
cudos por quilowatt da ponta. 

b) O funcionamento deverá ser fortemente automati- 
zado para evitar a necessidade de pessoal acompanhando 
permanentemente a exploração. 

c) Os encargos com a exploração serão muito atenuados 
se o funcionamento da central se fizer apenas durante cerca 
de 16 horas diárias (das 7 às 23 horas, por exemplo). Toda- 
via o funcionamento não permanente impede o uso do fri- 
gorifico doméstico eléctrico o qual, além de ser comodi- 
dade muito desejada pelos consumidores, contribui fortemen- 
te para aumentar a utilização da potência de ponta. 


Tendo em conta encargos que não se consideraram ao 
calcular o valor de €, (expressão 2) tais como 


— transporte do combustivel até à central; 
— pessoal atento ao funcionamento; 

— lubrificantes; 

— revisões; 

— reparações e substituição de peças; 


pode admitir-se que o custo da energia fornecida aos con- 
sumidores será cerca de 20 a 25 esc/k Wh. Este custo é hoje 
(1983) um dos padrões para avaliar o interesse de fontes 
alternativas. O seu valor, sem dúvida elevado, abre perspec- 
tivas a estas fontes de energia. 

À economia da produção de energia eléctrica pode, em 
princípio, ser melhorada se parte do calor dos gases de 
escape dos motores for recuperada, aquecendo água, de- 
pois utilizada para fins domésticos. 

Tendo em conta que os gases de escape estão a tempera- 
turas entre 200 e 300º € não se afigura dificil transferir boa 
parte da sua energia térmica para água que se acumula em 
recipiente termicamente isolado; o recipiente alimentará os 
consumidores domésticos. 

Esta recuperação merece ser cuidadamente estudada € 
demonstrada. 


3. ENERGIA DO VENTO 


O estudo dos aerogeradores está suficientemente avan- 
çado para ser possível uma primeira avaliação da energia 
que se pode extrair do vento [2]. [3], [4], [5]. 

De [2] extrai-se o Quadro 2 que indica qual a potência 
P, [kW] que se pode obter de turbinas cólicas com pás de 
diâmetro D[m], admitindo um factor de aproveitamento 
C, = 0,3, para diferentes velocidades do vento v [km .h!. 


QUADRO 2 — Potência de turbinas cólicas [kW] 


45 [km.h 1] 


D v= tt 20 30 
[m] 
2 0,012 0,1 0,33 11 
5 0,078 0,6 3 7 
IO 0,31 2.5 8,4 28 
20 1,24 9,9 33 [13 
40 5 40 134 452 
60 11,2 89 301 1017 
100 31 248 837 2825 


O Quadro 2 é obtido a partir da expressão 


Py=C,. v*. D*. 10, 35.109 [kW] (3) 


Como se sabe [3], [4] o factor de aproveitamento, 
sempre inferior ao limite de BETZ (C, limite = 0,593), ra- 
ramente ultrapassa o valor 0,45, encontrando-se dispositi- 
vos que não ultrapassam C, = 0,15. 

A expressão 3 mostra a fortíssima dependência da potên- 
cia Py com a velocidade v: Py varia com o cubo da velo- 
cidade. Desta dependência resulta que a energia cólica se 
adapta com dificuldade a abastecer consumidores que re- 
querem grande regularidade no fornecimento de energia. 
Estão nesta situação os consumidores domésticos que pre- 
tendem ver a televisão diariamente, precisam de iluminação 
todas as noites, têm o frigorífico permanentemente ligado à 
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rede. Por isso afigura-se que os aerogeradores, só por si, 
não podem assegurar o abastecimento das redes rurais 1s0- 
ladas. 

Diversas são as possibilidades para complementar a 
energia do vento: por exemplo armazenamento electro- 
químico ou grupos diesel-eléctricos. O uso simultâneo de 
aerogeradores e de células fotovoltaicas regulariza o dia- 
grama de produção mas não dispensa os equipamentos 
complementares que se referiram. 

A exploração conjunta destes equipamentos é pista a 
investigar e demonstrar. Todavia, desde já se retem, que a 
necessidade de complementar as fontes cólicas quando se 
destinam a abastecer redes rurais isoladas é condicionante a 
ter em conta ao escolher as interfaces entre o motor e a rede 
eléctrica. Situação muito diferente é a ligação de aerogera- 
dores a redes eléctricas de grande dimensão às quais estão 
ligadas outras fontes de grande potência, mas não é disso 
que se está a tratar. 

Uma grandeza que importa conhecer é a utilização que se 
pode esperar dos aerogeradores; alguns números começam 
a delinear-se baseados nas primeiras explorações à escala 
industrial. 

Para aerogeradores instalados em áreas em que a proba- 
bilidade acumulada P', [9%] de ocorrência de ventos com 
a velocidade v [m. seg!) ou maior, se aproxima da ex- 
pressão 


r v 
logio Po = 2-5 (4) 


e em que a caracteristica mecânica (P,, v) é a represen- 
tada na figura |, a utilização da potência nominal parece 
ser da ordem de 2000 horas por ano (a utilização de 2000 h 
é extraída de um trabalho elaborado pelo Assistente Esta- 
giário Eng. Fernando Nunes). Esta utilização só é possível 
se toda a energia cólica apropriável pelo aerogerador for 
utilizada. 

A apropriação da energia cólica com o rendimento má- 
ximo obriga a variar a velocidade angular consoante a ve- 


FIGURA | — (extraida de [2]) — 


curva | 


Probabilidade acumulada (P',, v). 
curva 2 Característica mecânica (Pr, v). 

Vo — velocidade de arranque. 

VI velocidade a que se define a potência nominal Pr im. 
Vmx velocidade de corte. 


locidade do vento [3], [2]. Regista-se, pois, que a energia 
cólica se adapta bem à obtenção de tensões alternadas, mas 
de frequência variável. Uma pista a investigar e a demons- 
trar é a redução na energia anualmente extraível dos aero- 
geradores quando são regulados a velocidade constante. 
Esta informação é importante para escolher a interface 
eléctrica a usar quando a energia cólica se destina a abaste- 
cer redes rurais isoladas. 

Contudo a irregularidade no fluir da energia cólica pouco 
afectará outros usos também importantes para os consu- 
midores rurais: bombagem de água, moagem de cereais, 
refrigeração e aquecimento e outras actividades importan- 
tes na agricultura ou que a complementam. Neste particu- 
lar o uso da energia cólica, embora sem passar pela forma 
eléctrica, tem secular tradição em Portugal. 


4. GERADOR FOTOVOLTAICO 


O estudo dos geradores fotovoltaicos está, também, sufi- 
cientemente avançado para se avaliar, com alguma segu- 
rança, a energia que anualmente poderão fornecer [7], 
[8], [9], [10]. 

Embora se façam grandes esforços para melhorar o ren- 
dimento das células fotovoltaicas, o interesse dos investiga- 
dores põe a tónica na redução do custo por unidade de 
potência — progresso que se afigura menos dificil. Con- 
tudo só grande redução nos custos actuais tornará Interes- 
sante o uso das células fotovoltaicas para alimentar redes 
rurais isoladas. 

De facto a utilização anual h, [horas . ano !] da po- 
tência de ponta P, das células fotovoltaicas é dada por 


h. = Poed X 8760:X 1 (5) 
em que P,eg é à potência média anual [kW] da radiação 


incidente e m o rendimento [p.u] do equipamento entre 
os terminais da célula e a rede eléctrica. A potência de 


ponta é definida para o feixe incidente convencional de 
1000 W/m* 

Para as potências médias anuais P, «4 (média das 8760 
horas) entre 150e 250 W/m? [8], que se obtêm nas zonas 
em que a conversão fotovoltaica tem mais interesse, a utili- 
zação anual da potência de ponta variará entre 920 h. 
ano”! e 1540h. ano ”!. Valores tão baixos tornam difícil 
a competição com outras fontes de energia. 

Para que o custo por quilowatt-hora não exceda | = 0,2 
SUSA/kWh (cerca de 26 escudos em Janeiro de 1984) o 
investimento por quilowatt de ponta 1, não deverá exceder 
lay < ha. f/1. Para encargos anuais de capital à = 0,2 
(isto é 20% por ano) obtém-se Loi < h, ouseja um inves- 
timento entre 920 e 1540 SUSA /kW,. 

Os preços actualmente praticados são cinco a dez vezes 
superiores. Todavia há razões para admitir que os preços 
dos conversores fotovoltaicos continuarão a descer rapi- 
damente [8]. 

A energia obtida nas células fotovoltaicas oferece vanta- 
gens sobre a obtida nos aerogeradores: as variações na 
energia afluente são muito mais lentas (nada ocorre compa- 
rável às rajadas do vento); a afluência da energia solar é 
muito mais regular e previsível; entre os paralelos 40º N e 
40º S a energia solar é uma fonte aproveitável, uniforme- 
mente distribuida, enquanto os ventos com velocidades 
médias aproveitáveis apenas ocorrem em pequenas áreas 
[4]; por isso, no futuro, poderá suportar programas na- 
cionais de fornecimento de energia eléctrica a áreas rurais 
isoladas, muito melhor do que a energia do vento. 


5. A QUALIDADE DA ENERGIA DAS FONTES 
ALTERNATIVAS 


Os aspectos da qualidade da energia eléctrica que são 
importantes para a alimentação de redes rurais isoladas 
são: 

a) garantia da continuidade do serviço; 

b) regulação da tensão; 

c) forma da onda; 

d) regulação da frequência. 


Outros aspectos, importantes para fontes ligadas a gran- 
des redes, tais como a estabilidade sincrona e as trocas de 
energia reactiva, não são relevantes para as redes rurais 
isoladas. Deve acrescentar-se que os equipamentos electro- 
domésticos fabricados em grandes séries, são construidos 
para alimentação em corrente alternada e por isso a rede 
rural deverá receber da interface uma tensão alternada, 
bem regulada, sinusoidal e com frequência constante. 

A garantia da continuidade do serviço é importante para 
os frigoríficos domésticos e para a iluminação. Menos im- 
portante para a televisão e para a rádio-recepção, admite 
interrupções não programadas, com duração de algumas 
horas, quando se trata dos outros electrodomésticos men- 
cionados na subalínea 1.2. Já se notou, a propósito da 
energia de origem cólica, que a garantia da continuidade do 
serviço pouco interesse tem para usos que são importantes 
nas comunidades rurais: rega, moagem e outros trabalhos 
complementares da agricultura e pecuária. Mesmo a refri- 
geração de produtos agricolas admite interrupções no for- 
necimento da energia eléctrica desde que o “frio” se arma- 
zene em massas de gelo ou por processos equivalentes. Isto 
que não é cómodo para usos domésticos, é realizável, com 
economia, nas indústrias agricolas. 


Aceita-se a regulação da tensão numa banda relativa- 
mente larga: variações de t+ 5% são, de um modo geral, 
aceitáveis; redução até — 10%, embora inconveniente, não é 


impeditiva do uso da maioria dos receptores. Regulações 
da tensão, dentro destas bandas, são fáceis de conseguir 
com as fontes de energia alternativas. 

Quanto à forma de onda a situação é diferente quando a 
energia passa pela forma mecânica, caso em que a obtenção 
de tensões sinusoidais com taxa de harmónicas mínima é 
solução natural, ou quando se usam células fotovoltaicas 
que nos seus terminais dão uma tensão continua. 

A frequência deverá manter-se numa faixa de + la 2%. 
Nas turbinas cólicas e também nas turbinas hidráulicas de 
pequena potência, a regulação a velocidade constante, com 
esta tolerância, não apresenta dificuldades técnicas, mas 
tem dois inconvenientes: a) o custo dos equipamentos de 
regulação; b) a redução da quantidade de energia que é 
fornecida pela turbina. 

Confrontando as exigências de qualidade com a quali- 
dade da energia fornecida pelos conversores primários 
(turbinas ou células fotovoltaicas) ressaltam algumas das 
orientações que devem presidir à escolha das interfaces e 
das suas características funcionais. 

Em primeiro lugar parece indispensável o recurso à acu- 
mulação da energia: são inevitáveis as calmarias que anu- 
lam qualquer possibilidade de extrair energia do vento du- 
rante 10 a 30% das horas do ano; as variações da potência 
da radiação solar não seguem o diagrama de carga; a irre- 
gularidade da afluência da energia hidrica é, em geral, 
muito grande, a tal ponto que só em casos particulares tem 
interesse o uso das centrais mini ou microhidroeléctricas 
para abastecer redes rurais isoladas, sem forte apoio diesel- 
eléctrico. 

Das acumulações de energia possíveis apenas se afigura 
com interesse para redes rurais isoladas a acumulação elec- 
troquímica. Adiante algo se acrescentará sobre este as- 
sunto; por agora retem-se que a interface deverá incluir 
acumuladores electroquímicos ou grupos diesel-eléctricos 
(que têm energia acumulada no depósito de combustivel) 
ou ambos os equipamentos. 

Do confronto da necessidade de tensões alternadas, com 
a tensão continua d- da pelas fontes fotovoltaicas resulta 
necessária a conversão continua-alternada o que sugere o 
uso de pontes conversoras. É sabido que as pontes conver- 
soras podem ser usadas como conversores da frequência 
variável em frequência constante. Se forem usadas com este 
fim podem evitar o uso de reguladores de velocidade nas 
turbinas, permitindo a apropriação das energias naturais 
com o rendimento máximo. Retem-se pois que, em qual- 
quer caso, as pontes conversoras são componentes a usar 
na interface. 

A tensão de saida das pontes onduladoras embora facil- 
mente regulável em amplitude e frequência, contem taxa de 
harmónicas incompatível com boa parte dos usos que inte- 
ressam aos consumidores rurais. 

Não é conhecida definição fundamentada das taxas de 
harmônicas que impedem o uso dos receptores que interes- 
sam às comunidades rurais. Parece evidente que as lâmpa- 
das incandescentes pouco ou nada são afectadas, que os 
receptores rádio podem ser muito afectados e que os moto- 
res eléctricos estarão numa situação intermédia. Todavia 
tudo isto é vago e requerem-se quantificações precisas. É 
outra pista que requer estudo e demonstração; pode, con- 
tudo, reter-se que as interfaces deverão incluir filtragens 
nas saídas com frequência estabilizada. 


6. INTERFACES 


Da alínea anterior pode concluir-se que uma interface 
genérica terá a constituição descrita pela figura 2. Esta in- 
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FIGURA 2 — Interface genérica, sem os órgãos de controlo. 


EÓL. 


CÉLULAS 
FOTOVOLTAICAS 


terface genérica permite analisar um modelo único, embora 
a potência e o modo de funcionamento dos diversos com- 
ponentes variem de caso para caso. Alguns desses compo- 
nentes (o grupo diesel-eléctrico, por exemplo) poderão até, 
ser suprimidos o que corresponderá a casos-limite da situa- 
ção modelada. 


Do gerador G, qualquer que seja, obtem-se, provavel- 
mente, frequência variável o que justifica a rectificação e 
subsequente ondulação, A irregularidade na afluência quer 
da energia cólica quer da radiação solar, justifica a acumu- 
lação electroquimica. A grande irregulariadade na afluên- 
cia da energia hidrica só pode ser compensada por grupos 
diesel-eléctricos. A pesada taxa de harmônicas na saída do 
ondulador e as trocas de energia reactiva justificam a filtra- 
gem e compensação €. 


Equipamentos complementares, variantes nas ligações e 
outros arranjos são, evidentemente, possiveis e a literatura 
especializada está cheia de sugestões quanto a interfaces 
possíveis. (Ver, por exemplo, o Capitulo 3 da referência 
[3]). Alguns dos arranjos que não se enquadram na tenta- 
tiva de sistematização representada na figura 2, têm mérito 
em situações particulares. Todavia afigura-se que as solu- 
ções redutiveis ao esquema genérico da figura 2, são as que, 
mais provavelmente se vão generalizar. 

Os componentes da interface descrita (gerador sincrono, 
gerador hipersincrono, rectificador, ondulador, acumula- 
dor electroquimico, condensador) são bem, ou pelo menos 
razoavelmente, conhecidos pelos electrotécnicos. O que é 
necessário investigar são as interacções nos sistemas em que 
coexistem e interactuam alguns ou todos esies componen- 
tes; e quais as características nominais óptimas para cada 
componente. Os principais condicionamentos resultarão do 
diagrama de carga e da qualidade de energia que se preten- 
dem na rede alimentada. 

Já se notou que do confronto entre a irregularidade da 
afluência da energia nas fontes primárias e as exigências de 
continuidade dos consumidores, resulta a necessidade de 
recorrer a acumulações por via electroquímica. Outros 
meios de acumulação de energia serão, certamente, possi- 
veis no futuro, mas em 1984, para redes rurais isoladas, a 
acumulação electroquimica ainda se afigura preferivel. 
Pode perguntar-se qual o tipo de acumuladores a usar. 

Os acumuladores “chumbo-ácido” usados nos veiculos 
rodoviários são construidos para satisfazer solicitações de 
grandes potências durante os curtísssmos tempos do arran- 
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que e fóra disso apenas funcionam como regularizadores 
(maus) da alimentação dos faróis, da ignição e pouco mais. 
Para estes fins bem específicos, desenvolveram-se determi- 
nados tipos de acumuladores de chumbo com electrólito 
ácido que se adaptam mal ao armazenamento e devolução 
da energia eléctrica nas redes rurais isoladas: por exemplo a 
energia armazenada por unidade de massa, a insensibili- 
dade à temperatura ambiente, a possibilidade de cargas e 
descargas rápidas são muito importantes nos veiculos e 
quase indiferentes nos usos estacionários; precisamente o 
contrário acontece quanto ao rendimento (energia) e ao 
custo inicial. 

Existem acumuladores electroquimicos que usam outras 
reacções e satisfazem melhor os requisitos das redes rurais; 
e estão a desenvolver-se outros tipos (por exemplo acumu- 
ladores secundários sólidos) que poderão vir a revelar-se 
ainda mais adequados. Todavia a assistência universal- 
mente estabelecida é a que serve os acumuladores “chum- 
bo-ácido” da tracção rodoviária; é imprudente contar com 
esta assistência para a inspecção, manutenção e reparação 
de outros tipos de acumuladores. 

É da experiência comum a importância que a assistência 
tem para a eficiência e duração dos acumuladores electro- 
químicos. Por isso justificam-se dúvidas, longe de estar es- 
clarecidas, quanto à linha de desenvolvimento que interessa 
seguir: ou aproveitar as novas soluções para a acumulação 
electroquimica e pode perguntar-se como se assegura a as- 
sistência em locais distantes dos grandes centros tecnológi- 
cos; ou tentar desenvolver tipos “chumbo-ácido” melhor 
adaptados às exigências das redes rurais. 


7. CRITÉRIOS DE VALORIZAÇÃO ECONÓMICA 


Na figura 3 a curva | representa o diagrama de carga 
sobre uma rede rural isolada, A curva 2 representa como 
varia no tempo a potência fornecida por uma fonte alterna- 
tiva. À potência pedida pela rede, P,, acima de P,, é for- 
necida por uma central diesel-eléctrica. 

A energia fornecida pela fonte é Wi = Wi + Wp. A 
parcela Wy carrega a bateria que restitui Wp = np. We 
sendo ng o rendimento global do armazenamento electro- 
químico. Ao consumidor é fornecida a energia Wj + Wp + 
+ Wp. sendo a parcela Wr obtida na central diesel-eléc- 
trica. Admita-se que no diagrama representado T = 24 horas 
e que o diagrama se repete todos os dias do ano. 


Se for aumentada a potência nominal da fonte alterna- 
tiva será menor a energia W, pedida à central diesel, mas 
será também menor a utilização da potência da fonte alter- 
nativa. Note-se que a capacidade da acumulação quimica 
cresce com a potência nominal da fonte alternativa. 

O grupo diesel-eléctrico terá a potência Ps > Ppnp;a 
potência Pr, justificar-se-á para garantir a segurança do 
abastecimento de cargas essenciais. Com a variação da po- 
tência nominal da fonte alternativa não varia a potência Pp, 
apenas varia a quantidade de energia Wp que a central 
fornece. 

Seja lpg| O investimento unitário [escudos / kW] no 
conjunto “fonte alternativa + bateria”, sendo upy a utili- 
zação do conjunto, referida à potência nominal da fonte, 


Wa + Wp 
Big E 56 ===. (6) 
F 
Se for 
W, 
tg =365=— (7) 
P.. 
F 
será 
Usg = Up [1 - KO -mng)] (8) 


Como se sabe up, mede o número de quilowatt-hora for- 
necido por ano à rede por cada quilowatt de potência no- 
minal da fonte. 

Se a parcela q. leg do investimento for paga inicial- 
mente (t = 0) e se os encargos do capital em divida forem 
anualmente y (| — a) leg. durante o tempo t,, os dis- 
pêndios durante o tempo t,, actualizados à taxa a e de- 
flacionados à taxa b, serão 


= ty vY(l-a) les 
ol a + EE 
di = 1; (| +a+b)! 


a que há a acrescentar os encargos com a exploração (revi- 
sões, reparações e substituições de componentes) 


T R, 


DR toe 


t; (Ll+a+b)! 


FIGURA 3 
— Diagramas, 


em que R, é o encargo de exploração (por unidade de po- 
tência P,.) pago ao fim de cada intervalo de tempo... Ter 
a vida física do equipamento. 

As receitas podem ser valorizadas pelo somatório 


li fo Upg (1 + e)! 


t=t “(+a+b! 


(1 +c)! 
= fo ue li=KOa)]——— (9) 
en TE E sas bi! 


em que c éa taxa de deriva dos preços da energia diesel fp; 
deverá acrescentar-se o valor final V; do equipamento, de- 
corrido o tempo 7, (também actualizado e deflacionado), 
[escudo . kW de P 


Vi 


(l+a+b)r 


Usando o critério do “Balanço Actualizado”, a utilização 
ur que anula a soma algébrica de receitas com despesas é 


| | tp o (l+a+b)! ty (L+a+b) 
(10) 
MW | | 
(U+a+b)” Tr tqi+o! 
fo Tiga) E ——— 
ni » ty (I+a+b)! 


O exame do diagrama de carga permite estabelecer a liga- 
ção entre a potência nominal P, da fonte e a sua utiliza- 
ção up (figura 4). 
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Fixado up pelo critério acima (expressão 10) fica definida 
qual a potência P,. a instalar. 

O valor a atribuir à energia fornecida pela fonte alternativa 
(valor de 1, nas expressões 9 e 10) poderá ser o valor dos 
encargos marginais da produção diesel. Isto corresponde à 
hipótese feita, das variações de P, não afectarem a potência 
do grupo diesel P,,, ou seja que em qualquer caso P',, < P,,. 

Obviamente do exame dos diagramas da figura 3, pode 
extrair-se a função uUpg (Pp). Nesse caso a expressão IO será 
alterada para dar imediatamente upp. 
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